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Abstract: Code optimization is an important area of research that has remarkable con- 
tributions in addressing the challenges of information technology. It has introduced a new 
trend in hardware as well as in software. Efforts that have been made in this context led to 
introduce a new foundation, both for compilers and processors. 

In this report we study different techniques used for sequential decomposition of mappings 
without using extra variables. We focus on finding and improving these techniques of compu- 
tations. Especially, we are interested in developing methods and efficient heuristic algorithms 
to find the decompositions and implementing these methods in particular cases. We want to 
implement these methods in a compiler with an aim of optimizing code in machine language. 
It is always possible to calculate an operation related to K registers by a sequence of assign- 
ments using only these K registers. We verified the results and introduced new methods. We 
described In Situ computation of linear mapping by a sequence of linear assignments over 
the set of integers Z and investigated bound for the algorithm. We introduced a method for 
the case of boolean bijective mappings via algebraic operations over polynomials in GF(2). 
We implemented these methods using Maple. 


Key-words: Linear algebra, number theory, boolean mappings, polynomials, sequential 
computation, compiler / processor / memory optimization, in-place algorithms, algorithm 
performance, multi-objective search. 
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Décomposition séquential des opérations pour 
l'optimisation des compilateurs 


Résumé : L'optimisation de code est un domaine de recherche en plein essor et dont 
les contributions pour faire face aux défis inhérents aux technologies de l'information sont 
considérables. Elle a introduit une nouvelle tendance dans le matériel ainsi que dans le 
logiciel. Les efforts qui ont été réalisés dans ce contexte ont conduit à introduire de nouveaux 
fondements, à la fois pour les compilateurs et les processeurs. 

Dans ce rapport, nous étudions différentes techniques de décomposition séquentielle in situ 
d'applications. Nous nous concentrons sur la recherche et l'amélioration de techniques de 
calcul et sur le développement d'heuristiques efficaces pour trouver ces décompositions, et 
mettons en ceuvre ces méthodes de calcul dans un compilateur afin d'optimiser du code en 
langage machine. Il est toujours possible de calculer une opération reliée à K registres par 
une série d'affectations à l'aide de ces seuls K registres. Nous avons vérifié les résultats 
et introduit de nouvelles méthodes. Nous avons décrit le calcul in situ de la cartographie 
linéaire par une séquence d'affectations linéaires sur l'ensemble des entiers Z et étudié la 
complexité de l'algorithme. Nous avons introduit une méthode pour le cas des applications 
via des opérations booléennes bijectives algébriques sur les polynómes en GF(2). Nous avons 
mis en ceuvre ces méthodes en utilisant Maple. 


Mots-clés : Algöbre linéaire, théorie des nombres, fonctions booléennes, polynómes, 
calcul séquentiel, optimisation de la mémoire / compilateur / processeur, algorithmes ön 
situ, performance algorithmique, recherche multi-objectifs. 
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Background 


This chapter consists of three main sections. In the first section, we describe the basic 
concepts of compiler and processor and the necessary background which is required to un- 
derstand the terminologies. "The second section explains the motivations and objectives. 
The last section describes a number of applications of compilers and processors. 


1.1 Introduction 


In this section, we describe some basic concepts about compilers, processors, registers op- 
timization, code optimization, techniques for compiler optimization, examples of compilers, 
types of optimizations and classification of microprocessors. 


1.1.1 Compilers 


A computer programmer typically writes softwares using high-level programming languages 
e.g. C++, Java etc. People can understand these languages but not directly the proces- 
sors. Assembly language is the language that the processor can understand (the machine 
language). Two different pieces of code in assembly language could be equivalent in what 
they do, but perform this work using a different sequence of steps. 

For example, when adding three numbers 1, 2and 3 together, there are different ways, the 
computer could execute this. One way would be to add 1 and 2 together and then add 
3 to that result (1 + 2) +3. Another way to add the three numbers would be to add 2 
and 3 together, and then add 1 to that result (2 + 3) +1. A compiler has many choices in 
which specific implementation of assembly language it will choose in making the translation 
from the high-level programming language [43]. A high level language must be converted to 
a language that a processor can understand, this task is performed by the compiler. The 
compilers were introduced in the early 1950's. 
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Figure 1.1: Complier interaction 


A compiler (more generally, translator) is a software application that converts the code 
written by the programmer into machine language [18]. For example, 


e gcc : converts C/C++ programs to assembly / machine code 
e f2c : converts Fortran programs to C programs 

e latex2html: converts Latex documents to HTML documents 
e Java : converts Java programs to JVM byte code 


ps2pdf: converts Post Script files to PDF files 


1.1.2 Compiler Optimization 


The main purpose of compiler optimization is to reduce the time taken by the program 
to be executed. The process focuses also to minimize the usage of memory storage and 
power consumption. Compiler optimization enables the program to be more efficient and 
ultimately the whole process helps to increase the speed for compilation. Usually three types 
of optimizations are considered: 


e Global optimization seeks to reorder the sequencing of a program in order to eliminate 
redundant computations (moving invariant operations outside loop bodies, colapsing 
loops, etc.). 


e Register optimization adjusts the allocation of machine registers to variables and in- 
termediate quantities in such a way as to minimize the number of a register has to be 
stored and later reloaded. 


CHAPTER 1. BACKGROUND 1 


e Local (time) optimization seeks to adapt the code to exploit particular features of the 
machine architecture and to remove local mishandling such as loading a register with 
a value that it already contains. 


There exist some code optimization problems that are considered to be NP-complete, or 
even undecidable [26, 201. 


1.1.3 Compiler Optimization Techniques 


There exist many compiler optimization techniques [48] but the work that has been done 
for the classification of these techniques can be improved. These methods can be divided 
into two main categories i.e. static and dynamic. These two approaches developed in 
parallel with a small interaction and optimization like parallelization has been obtained. 
Recently, the program parallelization of existing code came to forefront due to the fact that 
the major industrials focus to multi-core based architecture [47]. 


Static analysis of references to scalable data structures e.g. arrays can be very help- 
ful, for the purpose of thread level parallelization. It helps in proving sufficient conditions 
to observe the independence about the behaviour of array indexes, solving linear integer 
equations with constraints [30], special cases of array data flow (privatization) and outright 
pattern matching. 

Pointer analysis and shape analysis are also the important techniques but they are not very 
successful techniques [47]. 

Compiler optimization techniques operate on three levels: machine dependent, architecture 
dependent and independent. Instruction-level sensitivities of a compiler is described in the 
machine dependent. Architecture dependent includes parts of a program that relate to the 
general hardware implementation, but not to a specific machine. Architecture independent 
focus on that aspects of program formulation that do not depend on a particular computer 
system or even on a type of implementation (e.g. pipeline processing). Optimizations orig- 
inating in the academic and scientific community tend to be global, while, until recently, 
manufacturers have concentrated on local and machine-dependent techniques [53, 42]. 
Some Programming techniques take advantage of the optimizing compilers and the system 
architecture, e.g. BLAS, a library of Basic Linear Algebra Subroutines. The subroutines 
included in this library are able to provide significant enhancement in the performance of a 
program that is numerically intensive. 

ESSL, Engineering Scientific Subroutine Library. In fact, the ESSL library is an extension 
of BLAS library and includes high-performance mathematical routines for chemistry, engi- 
neering, and physics. 

PESSL (a Parallel ESSL) is another library that exists for SMP machines. Compiler options 
and the use of preprocessors like KAP and VAST, available from third-party vendors [13]. 
Several run-time techniques have been introduced by computer architects on the basis of 
hardware. This enabled the processor that it can run any ready instruction from an instruc- 
tion window [32, 35, 56]. Costly pipeline stalls has been prevented by bad execution and 
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detection of branches. Now the instruction can be stored in their dynamic execution order 
by using trace cache [21, 45]. Many run-time optimization has been introduced on the basis 
of hardware by latest microprocessors. 


1.1.4 Processors 


The processor, also known as the microprocessor (designed for microcomputers and micro 
controllers) or CPU (Central Processing Unit), is a complete computation engine that is 
integrated in a single chip. In fact, it is an integrated circuit, containing the arithmetic, 
logic, and control circuitry, and is used to interpret and execute instructions from a com- 
puter program. 

This integrated circuit, in combination with other integrated circuits that provide memory 
to store and execute the program, form a chip. 

Microprocessor registers used to hold temporary results, when the computation is being 
performed. Since these registers and microprocessors are made by the same technology, 
therefore there is no speed disparity between them. Moreover, these registers acts as pro- 
cessor memory. 

To improve the performance of microprocessor a small memory has been introduced between 
microprocessor and main memory. This small memory is called cache memory, it is expen- 
sive but fast and is first time introduced in IBM 360/85 computer. The first microprocessor 
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Figure 1.2: Microprocessor 
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was the Intel 4004, introduced in 1971. It was able to perform only subtraction and addition 
up to 4 bits at a time but everything was first time on a single chip. 
The microprocessors can be classified on the basis of: 


e the semiconductor technology of their design e.g. T'TL, CMOS or ECL. 
e the width of the data format (4, 8, 16, 32, 64-bit) they process 
e their instruction set e.g. CISC or RISC 


Due to the low power consumption, CMOS (complementary-metal-oxide semiconductor) 
technology is preferred to used in portable computers and in other devises that used batteries 
while TTL (transistor-transistor logic) is commonly used. 

When high performance is needed, the ECL(emitter-coupled logic) is used. For example, 
Older high-end mainframe computers, like the Enterprise System/9000 members of IBM's 
ESA /390 computer family, used ECL but current IBM mainframes use CMOS [6]. 
Four-bit devices are good for simple control applications and they are not so costly. 

CISC, (complex-instruction-set computer) processors, which have 70 to several hundred in- 
structions, are easier to program than RISC, (reduced-instruction-set computer) processors, 
but are slower and more expensive [2, 28, 52, 55]. 


1.1.5 Processor optimization 


How efficiently and effectively the processor executes instructions (provided in the form of 
a program designed by using some high level language) is determined by its internal design, 
also called its architecture. A processor can be considered as a combination of small blocks 
that organize to make a system. To help optimize the design space, a model is required 
that can predict the performance of a processor as a function of the delays of the underlying 
blocks. With such a model, one can evaluate how a change in the delay of a given module will 
affect the system's performance and can use this information to optimize a design. Because 
the design space is complicated, therefore, it may be difficult to know how changing the 
delay of a module will affect the overall performance of a processor [2]. 

Processor architect continue their efforts to improve the performance of processor every year. 
Some of the major techniques used by processor architects are the use of Wider data buses 
and registers, Floating point Units, Pipe lining and super scale architecture. 

As processor speed continues to increase faster than memory speed, optimization to use the 
memory hierarchy efficiently become ever more important. 

Blocking [24] or tiling [60] is à well-known technique that improves the data locality of 
numerical algorithms [1, 29, 31, 38, 40]. The improvement obtained from tiling can be far 
greater. Tiling can be used for different levels of memory hierarchy such as physical memory, 
caches and registers; multi-level tiling can be used to achieve locality in multiple levels of 
the memory hierarchy simultaneously [59]. The number of registers available on a processor 
and the operations that can be performed using those registers has a significant impact on 
the quality of code generated by optimizing compilers. 
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1.2 Applications 


In recent years, the incorporation of computers into a large number of devices e.g. laptop, 
palmtop, telephone etc., has been increased significantly. In designing such devices the 
important factors that must be considered are space, weight, power consumption. But this 
will lead to the limited amount of memory availability. 

A similar requirement is to use more and more sophisticated software in such devices, 
such as encryption software in telephones, or speech or image processing software in laptops 
and palm-tops. 

It is not always possible to run an application on such devices that require memory more 
than the available. This makes it desirable to try to reduce the size of applications where 
possible [49]. 

A fast computer program (as a result of compiler optimization) is not only useful for 
computer scientist and computer architecture, it affects the general public as well. Compiler 
optimization helps to increase the efficiency and capabilities of not only the sophisticated 
software but also increase the demand of newly introduced computer based devices. 
Ultimately, compiler optimization directly affect the computer based technology used in 
particular as well as in common life. 

For example, an improvement in computer programs for the medical community can affect 
all communities. An improved resolution of images obtained through scanning process has 
a direct affect on doctors and patients. An optimization will lead the life of the public to 
an ease. 

A number of soft wares are being used only for the purpose of research to observe the 
output of different results. A maple program can take more than a week to see the output 
of a program, but if the same program yield the result within one hour or less, imagine, 
how fast will be the conclusions based on these results. 


Microprocessors are the result of the semiconductor industrys ability to place an ever- 
greater number of transistors in a single integrated circuit. Optimized processors at the 
heart of mobile products enable communications, computing, and multimedia functions 
to be efficiently executed. Low power, high performance, and small form factor are key 
attributes. Open standards enable OS flexibility. 

Microprocessors also play supporting roles within larger computers as smart controllers for 
graphics displays, storage devices, and high-speed printers. 

However, the vast majority of microprocessors are used to control everything from consumer 
appliances to smart weapons. The microprocessor has made possible the inexpensive hand- 
held electronic calculator, the digital wristwatch, and the electronic game. 

Microprocessors are used to control consumer electronic devices, such as the programmable 
microwave oven and DVD player; to regulate gasoline consumption and anti lock brakes in 
auto mobiles; to monitor alarm systems; and to operate automatic tracking and targeting 
systems in aircraft, tanks, and missiles and to control radar arrays that track and identify 
aircraft, among other defence applications. Intel&)Core"" i7 processor is the most ever 
advanced desktop processor, introduced by Intel corporation recently in 2008. The Core 
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Diverse Market 


Figure 1.3: Processor Application 


i7 processor is the first member of a new family of Nehalem processor designs and is the 
most sophisticated ever built, with new technologies that boost performance on demand and 
maximize data throughput. The Core i7 processor speeds video editing, immersion games 
and other popular Internet and computer activities by up to 40 percent without increasing 
power consumption. 


Figure 1.4: Intel&)Core"" 17 processor 
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1.3 Motivations and objectives 


Since half a century, code optimization is considered to be an important area of research. 
It has introduced a new trend in hardware as well as in softwares. Efforts that have been 
made in this context led to introduce a new foundation, both for compilers and processors. 
New ideas have been introduced e.g. inter procedural whole program analysis, coloring- 
based register allocation, static single assignment form, array dependence analysis, pointer 
alias analysis, loop transformations, adaptive profile-directed optimizations, and dynamic 
compilation. 

Multi core processors have been introduced almost in all new computers. This multi core 
trend in the computer industry is forcing a new paradigm shift in compilers to face new 
challenges [50, 51]. The code of the program that is to be compiled has high effect on 
optimization. A large number of optimization has been introduced by modern compilers. 
These compilers have different effect on quality and size of code, time taken and energy 
consumption etc. [27]. 

In [7, 8, 9, 10, 11], it is shown that it is always possible to calculate an operation related 
to K registers by a sequence of assignments using only these K registers. Moreover, if this 
operation is linear or bijective, the number of assignments is at most 2k. In the general case, 
this number of assignments is at most 4k. 

For example, the bijection E on 3 bits defined by 


E(A, B, C) =(1+ A4 C + AB, A, B + AB + AC) 
can be calculated by a sequence of 4 assignments as given below: 
A:=A+B+C 
B:=A+B3+C 


Ci AO FAB 
A:=1+A+C+BC 


Such kinds of calculations generalize the traditional principle of the exchange of two bits 
A, B by the sequence: 


A:=A+B 
B:=A-B 
A:=A-B 


We are intended to study these methods and to improve these techniques of calculations, 
and especially to develop methods and efficient heuristic algorithms to find these decompo- 
sitions and implement these methods in compilers. 

We are interested in improving the bounds of these calculations, to obtain new methods 
of calculation for particular cases and implementing these methods in a compiler with an 
aim of optimizing code in machine language. 
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It is a fact that current computer architectures reach their theoretical limits of performance. 
However, it is still possible to gain performance of calculations by performing calculations 
in a new way. 

This implementation can be provided in hardware through the design of new processors, in 
the software, by optimizing compilers upstream (pre-level language compilers) and down- 
stream (post-machine language compilers). 

The current architectures use one or more processors equipped each one with a relatively 
low number of registers. 

However, these registers are constantly requested in the operations. İt is often advisable to 
minimize the use of registers. 
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| 
IN SITU Computations 


The conversion of input into output under well defined sequence of basic computational 
steps leads to the theory of sequential computations. A sequential program implements a 
mathematical function that maps a set of inputs to a set of outputs. These mathematical 
functions are well defined. The notion of computable functions has been introduced earlier 
by Church, Kleene and Turing. These functions are frequently used in untyped lambda 
calculus, recursive functions, and Turing machines [17, 33, 57]. 

Although, these basic models help in designing and reasoning for programming languages, 
domain theory and denotational semantics introduced by Scott and Strachey, and provide a 
global mathematical setting for sequential computation, building on top of the foundational 
theories [54]. This interconnects different programming languages and makes connection 
with the mathematical world of algebra, topology, and logic. It inspires the programming 
languages, type disciplines, and reasoning methods. 


2.1 In-place (In-situ) Algorithm: 


An in-place algorithm converts data structure using a minimal constant extra storage space. 
When such algorithms run, the input is overwritten by the output. For example, heap sort 
is an in-place sorting algorithm. 

An operation is said in-place operation if it does not alter the normal state of the system 
like a file backup can be stored over a running system without altering the speed of the 
system, while an in-place operation depends on the sophistication of the system. 

In order to improve cache performance, an algorithm or application should increase data 
reuse, decrease cache conflicts, and decrease cache pollution. A large amount of cache pollu- 
tion will increase the bandwidth requirement of the application, even though the application 
is not using more data [41]. 
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In the next example, two linear assignments are being reused and modified repeatedly to 
compute a linear mapping E. 


Example 1. The mapping E : (x,y) — (377x + 610y, 610x + 987y) 
can be computed by the following eight assignments: 


go —x — 2y 
y := —3T — y 
x := —rz—3y 
yı —3T — y 
zo =g 3y 
y := —3x — y 
x := x + 3y 

y := 2x — y 


Observe that, the sequence of eight assignments is actually the repetition of two opera- 
tions. The main advantage of such computation is the minimum usage of memory available. 


2.2 In Situ Decomposition of linear mappings: 


A linear mapping E : S” — S” can be decomposed into a sequence of mappings f : S" — 
S. The mapping E can be written as E : X :— AX, where A = la, il is a square matrix of 
size n, with a;,; € R and X is a vector with n variables. Each mapping f : S" — S could 
be considered as a linear assignment. These assignments can be computed by an in situ 
program that require no extra variables other than input variables. This in situ program 
can be simply written as follows: 
for i = 1 to k do 

Tp; = 04181 + Qi232 + +++ + Gin 
enddo 
where z4,2z5,::: , £n are the input variables and the output variables. The length of program 


is determined by the number 6 of assignments. 


Example 2. .Consider a mapping E that can be written as E : X := AX, such that 


21 2 3 5 8 
[xa _ [3 5 8 13 
X := ie and A := 5 8 13 21 
La 8 13 21 34 
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The given mapping E can be decomposed into the following sequence of 6 linear assignments: 


T1 := T1 T2 2x3 | 3x4 


Ta :— 41 T2 234 274 


£3 := 2x1 + 323 


£4 := £3 x1 + 272 
£1 := 33 — T1 — 272 


Ta := T3 — ZX] 


Observe that, the given mapping has been decomposed into a sequence of linear mappings 
in a way that, no extra variable other than the input variables is used. 


2.3 State of the Art:- 


In recent years, many computer scientists investigated the area of optimization including 
compiler optimization, processor optimization, mization, algorithm optimization etc. The 
optimization of specific linear algebra problems has been discussed on a large scale because 
such kinds of optimization effect on processor performance. Whaley and Dongarra discuss 
optimizing the widely used Basic Linear Algebra Subroutines (BLAS) in [58]. 

Chatterjee et al. discuss layout optimizations for a suite of dense matrix kernels in [16]. 
Park et al. discuss dynamic data remapping to improve cache performance for the DFT 
in [44]. Frigo, et al. in [22], which discusses the cache performance of cache oblivious 
algorithms for the matrix transpose, FFT, and sorting. Optimizing blocked algorithms has 
been extensively studied [36]. 


Most of the optimization algorithms depend on heuristics and approximations because 
several problems of optimizations are NP-complete. It may happen that a particular 
algorithm fails to produce better code or even worse. Code optimization transform a piece 
of code to make it more efficient without changing its output or side-effects [12]. 


An important optimization, affecting the performance of compiler code, is register 
allocation. Register allocation has been studied extensively in compilation and is a 
NP-complete problem. In 1981, Chaitin et al [14, 15] modelled the problem of assigning 
temporary variables to k machine registers as the problem of colouring, with k colours the 
interference graph associated to the variables. 


In general, register access is faster than memory access. Hence it is preferable to use 
register than memory whenever it is possible. When it is not possible to use register then 
some variables must be transferred to memory. This load/store operation has its cost. 
To avoid this cost some classical approaches have been introduced like in graph colouring 
algorithms [4]. 
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An iterated register coalescing algorithm, proposed by Appel and George [23] is a mod- 
ified version of previous developments by Chaitin et al. [14, 15], and Briggs et al. [5]. In 
these heuristics, spilling, coalescing (removing register-to-register moves), and colouring (as- 
signing variables to registers) are done in the same framework. These techniques are very 
useful in compiler optimization but still have to be revised to get better results. 

Burckel et al. [8] have introduced methods to design programs/electronic circuits, for 
performing any operation on k registers of any sizes in a processor, in such a way that one 
uses no other working memory (such as other registers or external memories). In this way, 
any operation is performed with at most 4k — 3 assignments of these registers, or 2k — 1 when 
the operation is linear or bijective. These methods are directly connected with processor 
and compiler optimization. 

In [7], it is proved that any mapping on 40,1)” can be computed by an in situ calculus 
without using extra variable. It is proved that three types of assignments are sufficient to 
perform these computations [9]. In [10], it is proved that any mapping E on (0, 1)” has a 
sequential computation in at most n? steps. In [11], it is proved that any linear mapping of 
n dimension can be computed with a sequence of 2n — 1 linear assignments. It is also proved 
that every mapping E on S" can be computed by an in situ program of length 5n — 4. The 
bound has been improved up to 4n — 3 when S is a power of 2. The maximal length of the 
program is 2n — 1 for bijective mappings [8]. 


2.4 Sequential Computation over R: 


In this section, we describe the proof that linear mappings on any field R admit sequential 
computation. Next, we explain, how the algorithm computes linear mappings sequentially 
over the field IR. To illustrate completely the ideas, we quoted different examples. We also 
quote some examples to explain, how the inverse mappings can be computed. We will also 
introduce an alternative approach to compute a linear mapping of dimension 2. We can also 
compute inverse mappings using this approach. We explain this approach by giving different 
examples. 


Proposition 1. Let E be a linear mapping on K”, where K is a field and n is a positive 
integer. There exists a sequence fi, fa,-** ,fn-1, fn, In-1,**" 592; gı of linear mappings 
from K” to K such that the program : 


21 := fi(z1,22,::- , fn) 
La :— fo(21,22,-- : Sta) 
Tn-1 :— fn-1(£1, a): qe ¿La 
Tn :— fe. (01,22, çə Un) 
Tn-1 := gain de 75 
T2 := go(1, To, + m 
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T1 I= gi(21, 22, ... ; Zn) 


performs the operation X :— E(X) for any X = (zx1,29,--- ,z4) in K”. 
In addition: 

For any i, gi will be x; or x; + v; for some j > i. 

E will be bijective => f;[i] 4 0 for any i. 

For any i, if fili] = 0 then gi = x; and fili] =0 for any j > i. 


Proof. The mapping E := E(X) can be considered as X := AX, where A is a matrix 
of coefficients i.e. aj; € K, for i,j = 1,2,---,n. Each row of the matrix AX represents 
the components of mapping E at current values of the variables (£1, £2, 9, za) and is 
considered to be a linear mapping . For example, the jth row ayı + aj2 + aya +... din, of 
matrix AX will be the E; component of mapping E. These linear mappings can be denoted 
by F; for i =1,2--- ,n. To manage the second sequence, the integers r; can be introduced. 
The first sequence f; can be computed by an iterative procedure as given below: 

For i from 1 to n, the following steps will be performed keeping r; = ? at the beginning. 


Case 1: If ai, = a Æ 0, then f; := F;, and x; will be modified to x; := ax; + A, where 
A does not depend on zi. For any aj; Z 0,7 > i, we will use the reference o ! (x; — A) 
to compute next mappings. 


Case 2: If ay, =a = 0, Vj > i, then f; :— F;, and x; := A, where A does not depend 
on zə. Since aj; = 0, Vj > 1, therefore no reference of x; will be required for the next 
steps. 


Case 3: If aj; = a = 0, and aj; Æ 0 for some j > i, then we can select aj; = 8 such 
that 8 4 0. We will perform the operation F; := F; — Fj i.e. subtract jth row from 
ith row so that 4 = —3 z 0, where & denotes the new value of a. In this way, we will 
be again in case 1 and by the hypothesis the value E; — Ej will be assigned to 2;. 
The operation F; :— F; — F} will introduce an assignment to the second sequence that 
can be obtained by adding Ej; to x; with r; :— j. 


Computation of the second sequence: 

The integers r;, that has been introduced, will be used to build the second sequence of 
assignments. 

Therefore from i = n to 1, the following steps will be performed iteratively. 


Step 1: If r; = i, then x; has been assigned in the first sequence and nothing to do 
more. 


Step 2: If r; = j such that j > i. Then because it results by the operation F; := F;—F;, 
therefore gj := x; + xj, but for à = n, this situation will not exist. 


The mapping E will not be injective in the Case-2 of the computation of first sequence but 
if aj; 4 0 for any i, then we can compute E”1 as well. 
Write the assignments, obtained, from last to first. 
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Compute x; := a^! (r; — A) from x; := ar; + A, (a 4 0). Use these references in the next 
assignments so that E”1 will be computed. 


2.4.1 Explanation: 


In this section, we will explain, the above proposition. We explain different cases by giving 
different examples for each case separately. 


Example 3. Consider a linear mapping E of the form 


E (21, £2, 23) (321 + 7x2 + 5x3, 8x1 + Avo + 913, 201 + xə + 623) 


We know that, for a vector X := (11, 12,23), the mapping E can be written as X := AX, 
where A denotes the matrix of coefficients. Therefore, for the given mapping E, matrix A 
will be of the form as given below. 


3 TS 

A=|8 4 9 

2.1 6 

Since a — 3 #0. Therefore, the first assignment is: 


11 :— 3414+ 722 + 523 (2.1) 
The initial value of x1 will be 
1 
3 (71 mi 7x2 = 513) 
Nou, for B = 8, we perform the operation 
T2 := 23 — Bar + B (o (xı — A)] 
on second row and for 8 = 2, we perform the operation 
T3 := 23 — Bar + B (oc (a, — A)) 
on third row of matriz A. 


1 
So, £2 := 821 + 4x2 + 9x3 — 821 + EC 7x9 zu 


8 44 13 


2 
and L3 := 2x1 + La + 6x3 224 + 3 (zı 7x2 523) 
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After performing this operation the matrix A will take the form 


3 7 5 
8 44 13 
A=|3 73 78 
2 11 8 
EE 


8 3 13 
a z z 
447 447 447 
The computation of third assignment: 
For, a= E and f = pe We perform the following operation: 


T3 := 23 — 23 + B (o! (v — A)) 


2 11 8 11 1118 3 13 
T3 := 371 3 Za + 373 | 3 T2 3 4471 442 1429) 
1 15 
z3 = irə + T” (2.3) 


Combining (2.1), (2.2) and (2.3), we get the required assignments, that are 


11 := 924 + To + 523 
8 44 13 


Ta :— 371 3 T2 3 T3 
me 
T3 := 472 4 T3 


Computing the Inverse of Matrix A: 


For a non-singular matrix A, we can also compute A~', by a sequence of assignments. We 
write, the sequence of assignments that can compute mapping E, from bottom to top, each 
assignment x; := ax; + À as x; := a7! (y; — A}. Applying this technique on the sequence 
of assignments 


zı 924 + Trg + 523 
8 44 13 


Ta := 371 3 Ta 3 T3 
1 Xx 15 
T3 := 472 4 T3 
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we get the corresponding sequence of assignments that can compute A^! and consequently, 
the mapping E”1. 


4 1 
£3 := — T3 — — 
Pec ng eb cud 
ə. 3,13, 
22705 a 
PE NEC: 
1: 3 1 372 373 
-1 37 43 
11 165 165 
. 2 8 13 
11 165 165 
İ 4 
15 15 


Example 4. Consider a linear mapping 
E (xy, £2, 23) — (5x2 + Taş, xo + 403, 311 + 223) 


For a vector X :— (xy, xə, xa), the mapping E can be written as X := AX, where A is a 
matriz of coefficients. Therefore matrix A will be. 


0 5 7 
A=|0 1 4 
3.0 2 


Since ayı = 0 =a. But aiz = 3 = 6. Therefore we have to subtract third row of matrix A 
from first row. The operation will be Ri := Ri — R3. 
After the completion of this operation, the matriz A will take the form: 


—3 5 5 
A=|0 1 4 
3 0 2 
Now the first assignment will be of the form: 
zı :— —311 + 5x2 + 5x3 (2.4) 


The initial value of xı will be 
1 
gcn + 513 + 523) 


Now, we will perform the operation 


Lo := Lo — Bay + Bia (rı — A)} 
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and 
T3 := 23 — Ba + B (o (21 — A)) 
The matriz A will take the form 


— 5 5 
A=|0 14 
—1 5 7 
The second assignment will be 
Lo := £2 + 4x3 (2.5) 


Now, the initial value of second assignment, xə will be 
£2 — 4x3 
And for the computation of third assignment: 
For, a = 1 and PB = 5, We perform the following operation: 
T3 := 23 — 23 + B (oc (v — A)) 
We will get the third assignment as 
£3 := —%1 + 523 — 1323 (2.6) 


The first assignment of second sequence that is obtained by the operation R4 := R4 — Ra will 
be. 


q1:— £1 43 
Hence the required assignments are 
Ly i= —321 + 5x2 + 5x3 
La := 22 + 4x3 


£3 := —%1 + 5135 — 1323 


Uy := Li T T3 
Computing Inverse of Matrix A: 


We can compute A~! by applying the same technique as described before i.e. by inverting 
the assignments and rewriting from bottom to top. 


T1 := Tı — T3 


1 
T3 :—= 13 (—21 + özə ET 23) 
_ 4 7 
T2 := 1371 132 
1 
T1 := 3 (=x1 + 5x2 + 523) 
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2 10 1 

39 39 

4 7 
Az — == 

13 13 

1 5 

13 13 


Example 5. Consider a linear mapping 


E (a1, 22, 23) — (3x1 + 522 + 723, 22 + 433, 273) 


For a vector X := (xy, xə, aa), the mapping E can be written as X := AX, where A is a 
matriz of coefficients. Therefore matrix A will be. 


3 5 7 
A=|0 1 4 
0 0 2 
Since ayı = 3 Z0 — a. Therefore the first assignment will be of the form: 
11 :— 3924 + 5To + T3 (2.7) 
The initial value of x1 will be 


1 
3 (71 Ll 5x9 EX 7x3) 


Since B = 0 both for second and third row, therefore, Matriz A will remain unchanged, and 
no reference of x1 is used to compute the second assignment. 
The second assignment will be 

La := £2 +423 (2.8) 


In computation of third assignment, we also do not need to use any reference of first or 
second assignment. Therefore, for such kinds of mappings the sequence of assignments can 
be written directly from the mapping. 


$3 :— 223 (2.9) 
Hence the required assigmments are 


q:— 311 + 5x» + 723 
Ta := T2 + 4x3 


13 := 2x3 


Computing Inverse of Matrix A: 
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The three assignments that can compute the inverse of matrix A can be written by applying 
the technique of inverting and rewriting the assignments from bottom to top: 


1 
T3 := 273 
T2 = xo — 4T3 
1 
zı = 3 (21 — 5x» m 723) 
1 5 13 
3 3 6 
A1-101 -2 
1 
00 3 


Example 6. Consider a linear mapping 
E (x1, T2, 23) =. (5x2 + 723, T2 + daş, 273) 


For a vector X :— (xy, xə, xa), the mapping E can be written as X := AX, where A is a 
matrix of coefficients. Therefore matrix A will be. 


0 5 7 
A=|0 1 4 
0.0 2 


Since ayı = 0 =a and B = 0 for all other cases, therefore we do not need to perform any 
operation, and the required assignments can be written directly from the given mapping or 
from the matriz A. 


£1 := 5x5 + 7x3 
La :— T2 + 4x3 


£3 := 223 
Computing Inverse of Matrix A: 


Since the matriz A is a singular, therefore A^! does not exist. It is also obvious that the 
assignment 
T1 i= 5x2 + 723 


is not an invertible assignment. 
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2.4.2 An alternative approach: 


Lemma 1. Every linear mapping E: (x, y) — (ma + ny, px + qy), where m, n, p, q € Z, 
can be computed by a sequence of at most 3 linear assignments with rational coefficients. 


Proof. Consider a sequence of three linear assignments: 


x := ax + by 
y :— cx + dy 
zo ez fy 


Where a, b, c, d, e, and f are rational In sequential computation, the linear assignments will 
take the form: 


x := ax + by 
y :— caz + (eb + d)y 
x := (ea + fca)x + (eb + fcb + fd)y 
The given Linear mapping E(2, y) (mx + ny, px + qy), can be viewed as a matrix of 


order 2 x 2 as below. 


(o 2) = Alsa) 


We can establish a system of four equations as given below: 


m = ea + fca 


n = eb 4- fcb + fd (5) 
pea 
g=bc+d 


Moreover, Determinant(A) =mx*q=—px*n=axexd 

— The product of three variables a, e and d should must be equal to the Determinant of 
matrix A. Therefore, a, e and d, will be the factors of Determinant of matrix A. 

The general solution of system (5) of equations is: 


“lı dE 5 qı O 
e e fp=m fp-m 

Bézout's identity states that, if two integers a and b are relatively prime then there exist 
u, v E€ Z such that au + bv = 1. 
Using Bézout's identity, we can construct a new matrix, such that we can compute the 
mapping. 
Suppose that m and p are relatively prime i.e. gcd(m,p) — 1. 
Using Bézout's identity, we can find u, v € Z such that um + vp = 1 


(2.10) 
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Multiply first row of the matrix A, by u and second by v, so that, u and v will be the 
common factor of entries of first and second row respectively, the new matrix will satisfy 
the required conditions. 

We will return to the original matrix by performing the operations. 


1 1 
y := (v and x := ə 
If m and p are not Co prime, then we can make m and p Co prime by extracting their gcd. 
We will solve the mapping as, for Co prime case and finally we will assign the GCD to a. 


'The above lemma can be explained by the following examples. 


Example 7. Consider the linear mapping 
E : (x, y) — (551 + 89y, 34x + 21y) 


such that matriz A will represent the coefficients. 


m n 55 89 
m iS 5 m E A 
Since GCD(55,89) = 1, Therefore, we can apply Bézout's identity to find integers u and v 
such that 55u + 80v = 1 


A number of solutions is possible for u and v. One of them is u = 13, v = —21 
Multiply first row of matrix A by u and second by v, we will get 


Aa (um un) _ 715 1157 
koup vq) \-714 —441 


IDEN pee IA a BE. ge ZE (2.11) 
e e Íp—m Íp—m 
we will have a = 1, 6 = 716, c = —714, d = 510783, € = 1, f = —1 


Hence the sequence of assignments is as given below: 


Using 


x := x + 716y 
1 
Un (-714x + 510783y) 
1 
2 


2.4.3 Computing the Inverse Mapping: 


We can also compute the inverse mapping using this alternate approach. We will explain 
the idea by giving different examples as follows. 
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Example 8. /n fact, the sequence of assignments that compute the mapping E represented 


by a matriz 
A ( 715 2 


—714 —441 
is given by 
x := x£ + 716y 
y := —714x + 510783y 
LILY 


The sequence of assignments that can compute the À! can also be obtained by inverting 
and rewriting the sequence of assignments used in computation of A. 


t= 
y := Ee {714r + y) 
510783 
x := x — T16y 
21 89 
fe 24323 39291 
34 55 


24323 39291 


We can compute the A^! by making some changes in the above sequence of assignments i.e., 
we will multiply the first column of A! by u and second by v. Therefore the sequence of 
assignments that can compute A^! is as under: 


T'=UXT+UXY 
1 
go x — 116y 


The A^! corresponding to matriz A is given below: 


21 89 
E 1871 1871 
AT! = 
34 55 
1871 1871 


LLC CC CC CCC CC CCC 


We give another example to explain the case, when the entries of the first column of matrix 
A are not Co prime 
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Example 9. Consider the mapping E: (x, y) — (25x + 13y, 35x + 21y), whose coeffi- 
cients can be expressed by a matrix. 


m mn 25 13 
4-( Jal 2) 


Since GC D(25,35) = 5, Extract ged from the first column, it will make the column co prime 


osa) 


Now GCD(5,7) = 1, Therefore Applying Bézout's identity 5u + Tv — 1, we can have a 
number of solutions, one of them is u = —4, v = 3 
Multiply first row by u and second by v, we will get 


pa (um un. —20 —52 
"vp vq) 121 63 


Solving the system of equations, we will have a= 1,b = 11,c = 21,d = —168,e = 1, f = —1 
such that the sequence of assignments for a = 1, is: 


x + lly 


8 
İl 


1 
A (21x — 168y) 
1 
—- 


Now just replace a = 5, The assignmenis for a — 5, are 
x := 5x + lly 


1 
yz (21x — 168y) 


r l(s-y 


ew 


2.5 Sequential Computation over Rings: 


In this section, we explain how a linear mapping can be computed by a sequence of linear 
assignments over rings specifically over the ring Z/NZ. We begin with a lemma of Em- 
manuel Thomé, next we will explain a proposition in detail. We will also explain the idea 
of computing a mapping over Z/NZ completely by giving different examples. 


Lemma 2. Let £1,--- , 24 be Co prime integers. Let N be an integer. There exists integers 
A2,*** , Àn such that zı +XE¡A¡2 € (Z/NZ)*, where (Z/NZ)* denotes the group of invertible 
elements. 
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Proof. Suppose that N is a prime power p". It is given that the integers z; are Co prime, 
therefore there exist an integer io (say) such that zio is Co prime to p. If zı is itself Co 
prime to p, then one can select ¿y = 1. If zı is not Co prime to p i.e. If xı is divisible 
by p, then zı + zio is Co prime to p. So the result holds when N = p" . For each prime 
power dividing N , we can therefore construct n integers AP ,--- , AP" such that A?” = 1 and 
Yi)? x; € (Z/P Z). 

Using the Chinese Remainder Theorem, these vectors combine into a global solution 
(1, A2, - , An) satisfying the required property. 


2.5.1 Assignment Matrices: 


Assignment matrices have been used to prove the result. An assignment matrix is actually 
the modified form of identity matrix having a row different from identity matrix. 


Definition 1. A matrix A is said to be an assignment matrix, if there exist an integer io 
such that for all row and column indices (i,j), one has either i = ig or Aij = öl. So, if A 
is a square matrix then A — I has at most one non-zero row, where I is the identity matrix. 


Example 10. Consider a square matriz A as given below: 


The given matrix can be decomposed into the following four assignment matrices under 
modulo 8 operation. 


Proposition 2. Let N be an integer. Any n xn matrix over (Z/NZ) can be written as the 
product of at most 2n — 1 assignment matrices. 


Proof. The result can be easily proved by the help of induction on the number of rows 
(representing by n) of the given matrix. Therefore if n = 1, the result is obvious. 

For n > 1, we can proceed as follows. Consider the first column of the given matrix À 
(say). Suppose that the azı is not an invertible element under modulo N. Let g represent 
the GCD of the first column i.e. g = (arı, dən, "nn ,@n1). We can construct an invertible 
element using the lemma as given above. We apply this lemma to make a combination 
of the coefficients of this column equal to g times an invertible element of (Z/NZ), with 
the constraint that this combination has its first multiplier equal to 1. This implies that 
the n x n matrix T defined by t;; = 6? for i > 1, and fr, = À; , where the multipliers 
1,12,---,), are obtained from the lemma above. Clearly the matrix T is an invertible 
assignment matrix, and the product TA has a coefficient at position (1, 1) which is equal 
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to g times an invertible element modulo N. Now assuming that as 1 € g(Z/NZ)* . Let 
n’ be the number of columns of A . Let G = diag(g,1,:-: ,1), and let A’ be the integer 
matrix AG! (A” has coefficients in Z because g is the GCD of the first column). We have 
a, € (Z/NZ*) . We form an n' x n' assignment matrix U defined by u;,; = 6) for i > 1, 
and uij = aj; . The matrix U is invertible modulo N (its determinant is a4 , ). The first 
row of the matrix A” = A” x UT! is equal to (1,0,--- ,0). Notice further that UG is an 
assignment matrix as well (even though not invertible modulo N ). Putting together the 
different results we have that A — T' x A" x (UG), where the matrix T may be omitted. 
Applying the result inductively on A" completes the proof. 


2.5.2 Explanation and Construction: 


For the given matrix A, if az 1 is invertible then matrix T will be the identity matrix of the 
same order as of matrix A and A” = AG !. 

Otherwise matrix T will be defined by t;; = ði; fort > landt,,; = A; 

i.e. the first row of the matrix T will consist of the multipliers 1,A1,--- , An, obtained by 
lemma given above. Matrix T may be a single matrix or T = T1, T5,:-- , Tn depending on 
number of prime factors of N. 

An invertible matrix L can be constructed by solving the system T'L — A. 

Le. Ly = T.A, Lə = TA,- , Ln = TrA, and A = (Ti*T5,--- , Ty) Ln. g is the GCD of first 
column of matrix L; obtained after T; transformation. So if A = (Ti x T5,--- , T;)L y. Then 
A'—L 7G, where Ly denote the final invertible matrix after 75,1 < i < n transformation. 
Matrix A will finally satisfy the relation A = (Ti x To,--- ,T,)A’UG. 


Example 11. Consider a mapping E : (11,12) — (2x1 + 312,511 +7x2). The coefficients 
can be represented by a square matriz 
2 3 
u- (5 7) 
Suppose that N — 8 — 25. 
For p" = 2, (2,2) 4 1, » m11 is not invertible. 


Since Az = 1, therefore 
. 1 1 
Matrix T= E i) 


6069-63 


a = —3 mod8 = 5,b = —4 mod 8 = 4,c = 5,d = 7 


New matrix 


Now, solving the system below 
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E dac AV LS: 0 25 4 5:5 
Mo E J il” İz z) mods = (1 3 


wem 9 )-( 3 
ve- 2 9-6) 


7 (10 11 PES 
TM ve- (5 7) mods = (5 3 


Hence the sequence linear assignments that can compute the given mapping is as follows: 


a := 5a + 4b 
b := a + 3b 
a:=a+b 


Example 12. Consider a mapping E : (11,12) — (2x1 + 312,521 + zə). The coefficients 
can be represented by a square matrix 


u-( 1) and N=6=2 x3 


mi, is not invertible with respect to p" = 2, (2,2) 21 
À = 1, and Az = 1, 


| (^1 A) (1 1 
2 ə. à 
Now, 
1 1) /a b) (2 3 
0 1)le d] 451 
=> a= —38mod6 = 3,b = 2,c = 5,d = 1 


New matrix 


Observe that lı 1 is not invertible, For p" = 3, (3,3) 4 1 
Ay = 1, and Az = 1, 
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Now 


=> a = —7 mod 6 = 5,b = 0,c = 5,d = 1 


New matrix 


T,T,M"UG — ES 1) mod 6 = ( 


5 1 


5 0 25 0 1 0 
( ) = Ge 1) modo = ( J 


Hence the sequence of linear assignments that computes the given mapping is: 


a := 5a 
b:=a+b 
a := a + 3b 


2.5.3 Computation of Inverse Mapping: 


In this section, we will explain, by examples, how to compute inverse of a square matrix, 
by a sequence of assignments over the ring Z/NZ. Consequently, the method will lead to 


compute E-!. 


Example 13. Consider a linear mapping 


E (xy, £2, xa) — (221 + 8x2 + 623, 321 + 1912 + 723, 524 +5%2 + 23) 
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We are interested to compute a sequence of linear assignment under modulo 9. 
Observe that, for a vector X := (11, 12,13), the mapping E can be written as X := AX, 
where A is a matriz of coefficients. Therefore matrix A will be. 


2 8 6 
A=|3 13 7 
5 5 1 
After performing modulo operation. 
2 8 6 
A=|3 4 7 
5 5 1 


The sequence of linear assignments that can compute mapping E, under modulo 9 operation, 
is given by 


11 := 2x1 + 8x2 + 623 


La := 6%, + L2 + 7T3 


£3 := Tri + 322 + T3 
Computing Inverse of Matrix A: 


We can get a sequence of assignments that can compute E^! and ultimately A7*. We obtain 
this sequence of assignments by inverting and rewriting the assignments (from bottom to top) 
that compute the mapping E. 


£3 := £3 — Tr) — 323 
Lo := Lo — 6x1 — 7x3 
T1 := 2-1 (zı = 8x2 = 6x3) 
In the last assignment, we will replace 27! by 5 because 2x5 = 1mod9, so that the last 
assignments becomes invertible modulo 9. 
Hence the sequence of assignments that can compute the inverse mapping under modulo 9 
is as given below: 
£3 := £3 — Tri — 323 
La := Tag — 624 = 723 


21 := 5 (£1 — 8x2 — 6x3) 


Consequently, we get the inverse of matrix A. 
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Verifying the product of assignment matrices: 


The matrix 
2 8 6 
A=|3 13 7 
5 5 1 


is decomposed into three assignment matrices under modulo9 operation, these matrices are 
as given below: 


2 8 6 1 0 0 1 0 0 
A=ļ|0 1 0] 42=|3 13 7|43=[0 1 0 
0 0 1 0 1 5 9 1 


0 
Clearly the product of matrices Aş x Az x A) = A. 


2.5.4 Inverses mod k: 


An integer a has an inverse mod k if and only if GC D(a, k) — 1. İt is not always true that 
a with a #0 (mod k) has an inverse mod k, e.g. 2% 0 (mod 4). Since 2 x 2 = 4 = 0 mod 
4. Thus 2 has no inverse. Otherwise, we could multiply both sides of 


2 x 2 = 0 mod 4 


by the inverse of 2 and get the false result 2 = 0 mod 4. If p is a prime number, then for 
each a % 0 mod p has a multiplicative inverse mod p. 


Existence of Inverse of Matrix A: 
Consider again the mapping 
E (xy, £2, xa) — (211 + 822 + 623, 321 + 1915 + 723, 524 + 523 + 23) 


as given in Example 2. This mapping can also be computed by the following sequence of 
assignments under modulo 12. 

zı 521 + 525 + 923 

xa := 3121 + 10x2 + 4x3 


£3 := 101 +42%3 


T4 :— 11 + 323 


By inverting and rewriting this sequence of assignments, we get the following sequence of 
assignments. 


Tı : = T1 — 319 
23 := 4 l (x3 — 21) 
x2 := 1071 (a2 — 324 — 473) 


xı := 5 l(x-—523— 9x3) 
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Observe that, the second and third assignments are not invertible because the integers 4 and 
10 have not their inverses under modulo 12 operation. Moreover, the Determinant of matrix 
A 
2 8 6 
LAL — 33 13 7 = 248 
5.5 1 


is not invertible modulo 12. Thus A^! modulo 12 is not computable by the sequence of 


assignments as given above. Hence, to compute inverse mapping the linear assignments 
should must be invertible. 
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In Situ Computation of Linear 
Mappings over Z: 


In this chapter, we describe In Situ computation of linear mapping by a sequence of linear 
assignments over the set of integers Z. We start in proving a proposition and then we explain 
the algorithm completely by giving different examples. 


Proposition 3. Let E be a linear mapping on Z”, where Z is the set of integers and n 
is a positive integer. There exists a finite sequence fi, fa, +--+: fn; Gn—1, *** > 92, 91 of linear 
mappings from Z" to Z such that the program : 


T1 :— filti, to, y Bix) 
T2 := (ai, to, ba) 
Tn i= Tat du : itn) 
Tn-1 := 9n—1(%1, La," əə T 
T2 := go(z1,223,--- , En) 
21 := g1(%1,%2,*** , Ln) 


performs the operation X := E(X) for any X = (11,%2,*** ,tn) in Z^. The sequence 
Qn: 9n—1: 77^. 92, 91 contains exactly n number of assignments. 


Proof. The mapping E := E(X) can be considered as X :— AX, where A is a matrix 
of coefficients i.e. aj; € Z for i,j = 1,2,--- ,n. Each row of the matrix AX represents 


a component of the mapping E at current values of the variables (x1, 22,--- r4) and is 
considered to be a linear mapping . For example, the jth row ayı + aj2 + aya +... din, of 
matrix AX will be the E; component of mapping F. 


Case 1: If aj; < 0, for j > i, then multiply jth column by —1, so that aşı > 0 and x; 
will be modified to x; :— —2;. Similarly, if a;; < 0, then ith column will be multiplied 
by —1, so that as; > 0 and x; will be modified to x; :— —x;. 


Case 2: If aj; > aij, forj > i, then the program will perform the operation C; :— 
C; — Cj, and x; will be modified to x; := ax; + A, where A does not depend on zi. 
C; denotes the ith column of matrix A. 


Case 3: If aij > au, forj > i, then the program will perform the operation C; := 
C; — Ci, and x; will be modified to x; :— ax; + A, where A does not depend on gp. 


Case 4: If aj; = 0, then the program will perform the operation C; :— C; — Cj, and x; 
will be modified to x; := ax; + A, where A does not depend on zi). 


The whole procedure works iteratively until the matrix A modified to be a triangular matrix. 
After the completion of this procedure the matrix A will become a lower triangular matrix 
L (say) as given below. 


lii 0 0 0 0 0 
loi ləə 0 ə 0 0 0 
006 7 0 0 
L = . . . . ” 
In—2,1 In-22 In—2,3 ... İn—2,n—2 0 0 
En In-112 Ín-13 +.» ln-1m-2 In-1yn-1 0 
la ln2 las ss ln n-2 lani İnn 


The last n assignments can be obtained directly from the matrix L and that will be of the 
form. 


Tn = Ini ®1 | ln212 5323 mu Inn 2Tpn-2 + lnn 17n—1 + lanTn 
Xn—1:— ln—1,171 HP In—1,2%2 a In—1,303 Pl 1,n-22n—2 + In 1n-1Un-1 
Tn-2 i= ln-2171 EE lj, 2,272 “İr In—2,3%3 e In—2,n—2€n—2 


23 := 13141 + 13222 + 13303 


To := lg121 + l22£2 


Xa := l11%1 
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Let A and V denote the modified form of the matrix A and vector V after performing any 
operation. After each operation it should must preserve and compute the original mapping 
effectively i.e. 
MV = MV 

If T' denotes the transformation matrix then 

M = MT 

Very 

=> MV = M(TT-V = MV 
Observe that, to compute mapping E the last n assignments are not necessarily to be 
invertible, because in this last sequence of assignments, the reference of initial assignments 
are not required to be used for the next assignments. 
The similar linear assignments can be combined to a single linear assignment that have the 
same effect. 


In general, if the assignment x; :— ax; + A repeated k times then we can write it as a single 
assignment of the form 


Ti = afri + A (o7! E of? +... + a? od o) 


We can modify the algorithm to avoid from repeated similar assignments. 


3.0.5 Explanation and Construction: 


Consider a linear mapping 


E (11, £2, £3) (5x1 — 3x2 + 523, 321 — Taş, 401 + 8x2 + 1323) 


We are interested to compute the given mapping by a sequence of linear assignments such 
that the coefficients of these assignments are integers. The coefficients of the given mapping 
E can be taken as a square matrix A: 


D —3 5 
Ai=|3 0 -—7 
4 8 13 


We shall compute the mapping E by applying the procedure, keeping in mind different cases, 
as described above in the proposition 2. Since a4» < 0, therefore, we perform the operation 


C2 :— —C2 because we are in the first case. The linear assignment corresponding to this 
operation is zə :— —zə and the matrix A will be modified to the form 
5 3 5 
A:—|[3 0 -7 
4 —8 13 
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Since a11 = 5 > a12 = 3, therefore, we perform the operation C1 :— C1 — C2 because we are 
in the second case. The linear assignment corresponding to this operation is 2 :— 11 + £2 
and the matrix A will be modified. 


2 3 Ə 
Ai=|3 0 -7 
12 -8 13 


Now, a12 = 3 > ay, = 2 therefore, we are in the third case and we perform the operation 
C2 := C2 — Ci. The linear assignment corresponding to this operation is x, := 11 + 12 and 
the matrix A will be modified. 


2 1 5 
4:-13 -3 —7 
12 —20 13 
Since a1, = 2 > qız = 1, therefore, we perform the operation C, :— C, — C2. The linear 
assignment corresponding to this operation is 13 := gi + x2 and the matrix A will be 
modified. 
1 1 5 
4:-16 -3 -7 
32 —20 13 
Since dış = 1 > aşı = 1, therefore, we perform the operation Cy := C2 — C1. The linear 
assignment corresponding to this operation is 11 := x, + x2 and the matrix A will be 
modified. 
1 0 5 
A:=|6 -9 -7 
32 —52 13 


Observe that a19 became zero, We will continue this process until dış = 0. 
Since dış = 5 > ayı = 1, we perform the operation C3 := C3 — C4. The linear assignment 
corresponding to this operation is x, := zi + 23. The matrix A will be modified. 


1 0 4 
Am i6 -—9 -13 
32 —52 —19 


Since dış = 4 > ayı = 1, we perform the operation C3 :— C3 — C1. The linear assignment 
corresponding to this operation is x, := 11 + z3. The matrix A will be modified. 


1 0 3 
4:-16 -9 -19 
32 —52 —51 


36 


CHAPTER 3. IN SITU COMPUTATION OF LINEAR MAPPINGS OVER Z: 3 


Since dış = 3 > ayı = 1, we perform the operation C3 :— C3 — C1. The linear assignment 
corresponding to this operation is x, := zi + 23. The matrix A will be modified. 


1 0 2 
Am 16 -9 —25 
32 —52 —83 


Since dış = 2 > ayı = 1, we perform the operation C3 :— C3 — C4. The linear assignment 
corresponding to this operation is x, := zi + z3. The matrix A will be modified. 


1 0 1 
A:=|6 -9 —31 
32 —52 —115 


Since dış = 1 > aj; = 1, we perform the operation C3 :— C3 — C4. The linear assignment 
corresponding to this operation is x, := 11 + z3. The matrix A will be modified. 


1 0 0 
A:=|6 -9  -37 
32 —52 —147 


Now the first column will remain unchanged throughout the next operations, we will focus 
on second diagonal entry. 


Since a23 = —37 < 0, we perform the operation C3 := —C3. The linear assignment corre- 
sponding to this operation is 13 :— —a3. The matrix A will be modified. 
1 0 0 
A:=|6 -9 37 
32 —52 147 
Since a22 = —9 < 0, we perform the operation C2 := —C2. The linear assignment corre- 
sponding to this operation is x2 :— —x2. The matrix A will be modified. 
1 0 0 
A:= | 6 9 37 
32 52 147 


Since a23 = 37 > a22 = 9, we perform the operation C3 := C3 — C2. The linear assignment 
corresponding to this operation is 12 := 12 + 23. The matrix A will be modified. 


1 0 0 
A:=|6 9 28 
32 52 95 


since a23 = 28 > a5» = 9, we perform the operation C3 :— C3 — C2. The linear assignment 
corresponding to this operation is 12 := zə + zs. The matrix A will be modified. 


1 0 0 
A:=|6 9 19 
32 52 43 
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Since a23 = 19 > a22 = 9, we perform the operation C3 := C3 — C2. The linear assignment 
corresponding to this operation is 12 := 12 + z3. The matrix A will be modified. 


1 0 0 
A:=|6 9 10 
32 52 —9 


Since a23 = 10 > a22 = 9, we perform the operation C3 := C3 — C2. The linear assignment 
corresponding to this operation is 12 := zə + z3. The matrix A will be modified. 


1 0 0 
A:=|6 9 1 
32 52 —61 


Since a22 = 9 > a23 = 1, we perform the operation C2 :— C2 — C3. The linear assignment 
corresponding to this operation is x3 := zə + 13. The matrix A will be modified. 


1 0 0 
A:=| 6 8 1 
32 113 —61 


Since a22 = 8 > a23 = 1, we perform the operation C2 := C2 — C3. The linear assignment 
corresponding to this operation is 13 := £2 + 23. The matrix A will be modified. 


1 0 0 
A:=| 6 7 1 
32 174 —61 


Since a22 = 7 > a23 = 1, we perform the operation C2 := C2 — C3. The linear assignment 
corresponding to this operation is 13 := 12 + 23. The matrix A will be modified. 


1 0 0 
A:=| 6 6 1 
32 235 —61 


Since a22 = 6 > a23 = 1, we perform the operation C2 := C2 — C3. The linear assignment 
corresponding to this operation is 13 := 12 + 23. The matrix A will be modified. 


1 0 0 
A:= 6 5 1 
32 296 —61 


Since a22 = 5 > a23 = 1, we perform the operation C2 :— C2 — C3. The linear assignment 
corresponding to this operation is 13 := 12 + 23. The matrix A will be modified. 


1 0 0 
Ami6 4 1 
32 357 —61 
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Since ag2 = 4 > a23 = 1, we perform the operation C2 := C2 — C3. The linear assignment 
corresponding to this operation is 13 := 12 + 23. The matrix A will be modified. 


1 0 0 
A:i=|6 3 1 
32 418 -—61 


Since a22 = 3 > 493 = 1, we perform the operation C2 := C2 — C3. The linear assignment 
corresponding to this operation is 13 := 12 + z3. The matrix A will be modified. 


1 0 0 
A:=|6 2 1 
32 479 —61 


Since a22 = 2 > a23 = 1, we perform the operation C2 := C2 — C3. The linear assignment 
corresponding to this operation is 13 := 12 + 23. The matrix A will be modified. 


1 0 0 
A:= | 6 1 1 
32 540 —61 


Since a23 = 1 > a22 = 1, we perform the operation C3 :— C3 — C2. The linear assignment 
corresponding to this operation is 12 := 12 + z3. The matrix A will be modified. 


1 0 0 
A:=|6 1 0 
32 540 —601 


The matrix A converted into a lower triangular matrix. The last three assignments 
will be written directly from this lower triangular matrix. These three assignments are: 
13 := 322, + 540z2 — 601273, £2 :— 67, + x and zı := 21 

Hence the given mapping E can be computed by the sequence of assignments as given below: 


T1 = T1 T T2 
Ta = Tı T T2 
T1 = T1 T T2 
Tı = Li T T3 
Li = Li T T3 
Li = Li T T3 
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T2 = T2 + T3 
23 = 32x, + 540%2 — 60123 
Ta = 621 + 2 


T1 = 21 


3.0.6 Computing Inverse mapping over Z: 


We can compute the inverse mapping E”1 over Z by inverting each assignment in the 
sequence (that can compute mapping E) and rewriting from bottom to top. We explain by 
the following example 


Example 14. Suppose that E be a linear mapping and is defined as 
E (x1, 22) — (8x1 + 1325, 1321 + 2123) 


The matriz of coefficients A is 


The given mapping E can be computed by the following sequence of assignments 
Xi := —21 — 222 


T2 := 3121 +22 


11 := 11 — 3922 


T2 :=2x1 + 22 


Now, to compute E~!, invert each assignment and rewrite from bottom to top. We get the 
following sequence of linear assignments that can compute ET. 


£2 :— £2 — 2x1 
“un 3x2 
Xo: £2 — 9x1 
21 := —X1 — 272 


Consequently, we can compute A7!, i.e. 


he Voir 13 
a “7 5 
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Existence of Inverse matrix 47): 


The computation of inverse mapping is not always possible over Z by a sequence of linear 
assignments. Consider the following example 


Example 15. Consider a mapping E 


E (a1, xə, v3) (221 + 312 + 5x3, 321 + 4x2 — 723, 8x2 + 1323) 


The mapping E can be computed by the following sequence of assignments. 


21 = T1 T T2 
zı = xı + 973 
T3 = —X3 
£2 = T2 + 773 
£3 = 24 x x) — 16 x zə + 219 * 23 
T2 = T2 
T1 = ZX] 


The matrix of coefficients A for the mapping E can be written as 


2 3 5 
4A4:-13 4 -7 
0 8 13 
and 
108 1 -4 
A71:— — | —39 26 29 
2191 24 —16 -1 


1 
But 219 € Z. Observe that the assignment 


v3 = 24 x xi — 16 x zə + 219 * 23 


is not invertible. Therefore the inverse mapping E-! can be computed by inverting and 
rewriting the above sequence of assignments. Moreover E-1 is computable by the sequence 
of assignments iff E is a bijective mapping. Hence the given mapping E is not a bijective 
mapping. 
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3.0.7 Investigating Bounds for the number of assignments: 


We are interested to find the minimum number of assignments required to compute a given 
linear mapping over Z. We proceed by developing different relations between mappings that 
can help to find the minimum number of assignments and we investigate through finding 
different counter examples. We are giving a proof below that shows that six assignments 
are not sufficient to compute any linear mapping over Z?. 


Lemma 3. İt is not possible that every linear mapping 
E: (a, y) — (mz + ny, pz + qy) 


, where m, n, p, q € Z, can be computed by a sequence of at most 6 linear assignments 


x := az + by 
y :— cx + dy 
xz :=exr + fy 
y := gx + hy 
x := ix + jy 
yo kz + ly 


where a, b, c, d, e, f, g, h, i, j, k, andl € Z, 
Proof. Consider a linear mapping E(x, y) — (461x + 286y,353x + 219y), that can be 


viewed as a matrix of order 2 x 2 as given below. 


461 286 
G 0 m. 


We can establish a system of four equations(by evaluating assignments) as given below: 


461 =a(ie+ife+ jge + jgf c + jhe) 

286 = ieb+ifcb+ifd+ jgeb + jgf cb + jg fd + jhcb + jhd 

353 = 461k + [gea + lg f ca + lhca (2) 
219 = 286k + lgeb + lg f cb + lg fd + lhcb + lhd 


Determinant(A) = mq — pn = aedihl 

=> The product of six variables a, d, e, h, 1, and £ should must be equal to the determinant 
of matrix A. 

Since, the determinant of matrix A is equal to 1. Therefore, a x e* də h: il — 1 and each 
of six variables can take only the value (1, —1). 

There are P(6,2) + 2 = 32, cases of assigning values to a, d, e, h, i, and l. 

Let us consider one case that isa=1,e=1,d=1,h=1,i=1,l=1. 


43 


We will have equations: 


460 = fe+ jg + 3gf c - jc 
286 = 461b + f + jgf +i 
353 = 461k + g+gfc+c 
218 = 286k + 353b — 4610k + gf 


There are three possible general solutions for system (3) of equations: 


The first possible solution is of the form: 


461 1 286 
B a IUS b — 
gt) jif =f, 4617 7 461 
hod —460j + 461f 13539? 1 
— 461 j? "dec mu 


'The second possible solution is of the form: 


460 " 1 —460 + 353j 


juge uec ; 
1 .. 286 
b= ) + 3 = J,C— 
f-t 4677 770 


The third possible solution is of the form: 
1 461g4-c— 353 — 353/g 


J—59—9 461 14 jg 
1 jg — 460 + 286c —460 + jg + jc 
c= c,b= v = . 
461 c c(1 + jg) 


Consider the first possible solution: 


1 
Suppose that g=--€Z = j=10or —1 
J 
Now, If 7 = 1 then 


g=-1€Z, but b= | 


Similarly, If ; 2 —1 then 
1 2386 287 
=1€Z, butb= = Z 
2. .270020: 


Therefore integer solution is not possible in the first solution. 
In the second possible solution: 


460 , 1 —460 + 3534 
j 461 j 
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—g + 353 


k b itt k= 
can be written as 461 


Observe that if k € Z then —g + 353 > 461 
=> g € 353 — 461 => g € —108 


But, since g is a divisor of 460, and, g can take only the value (—115, —230, —460) because 
g € —108. But for any of these values (—115, —230, —460) k £ Z. 
Notice that, 


230--353 6883 

if g 230 then k 461 461” 
1154-353 4468 

if g = —115 then k = = Z 

ig — 461 461 * 
460 +353 813 

if g 460 then k 461 181 $ 


Therefore integer solution is not possible using second possible solution. 
Now, in the third possible solution: 


1 jg — 460 + 286c —460 + jg + jc 
and f — : 
461 c c(1 + jg) 


Since c divides both jg — 460 + 286c and —460 + jg + jc, therefore, if b, f € Z then there 
should exist kı, kə € Z such that 


b= 


460 — jg — jc = cki (3.1) 
and 
460 — jg = c(286 — 461k2) (3.2) 


By (3.1) and (3.2) 
c(j + ki) = c(286 — 461k2) 


— j + kı = 286 — 461k, 


—286 + kı +j 
=> k —o——  ””” ””. 
? 461 
If ko € Z then —286 + kı + j > 461 
— kı +j > TAT (3.3) 
460 — ; 
But form equation (3.1) c = 27: 
jth 
and if c € Z then 
j + kı € 460 — jg (3.4) 
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By (3.3) and (3.4), we can write 
460 — jg = 747 => jg = —287 


=> For j and g to be integers, they should must be factors of 287 and 


(+1, +7, +41, +287} 


are the possible factors of 287. 
Let, for example, g — —287 then ; — 1 


ki T46 kə —1 C 1 


—460 — 287 + 1 34 
1 — 287 13 EZ 


But, then, f = 


For all divisors of 287, f £ Z. 
The other 31 possibilities of assigning values (1, —1) to a, d, e, h, i, and l also do not work 
and can be proved by the similar arguments. 


We have established a useful relation between different mappings with coefficients taken 
from Fibonacci sequence. Next, we explain a proof related to these relations. 


Definition 2. We will define a Fibonacci-like sequence as 
En = Fn-1 + Fn-2 
where Fo = F) = 1. 


Definition 3. We will define a Fibonancci-like matrix to be a matrix in the form 


Pn F, 
Fy Fuga 


o PU Faa ER 
E MEO Fu 
n F3 k Fin+1 P4n42 = Fanta Pings 
Py Fa Fan+2 Pang Fan+s Fang 
Lemma 4. Let E, : (x, y) — (Fa-ız + Fay, Fax + Fryiy) be the mapping on Z2, where 


F is the Fibonacci number. The mapping Eaxy2 can be computed with 2k + 2 number of 
assignments, where k = 0, 1, 2,--- n. 


We define also the relations 
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Proof. For n= 2, Es = (Fix + Foy, Fox + F3y) can be expressed as 
1 2 
#-( 3) 


The assignments for Aş = { 280 


y := 2z — y 


Li=xX+2y 
The assignments for Aş = 4 y:= 3x — y 
y:;—o—rc-y 
Observe that 
2-8 1 2 8 13 
As = (5 JE 3) =40= (i i 
gz := x + 2y ən 
y := 3x — y ao 
The assignments for Ag = 4 t:=-+Y = y= y 
v= 2 + 3y 
x := x + 2y Ee 
y := 27 — y 20 y 


M 2 3 8 13 55 89 
Similarly, Aş — E I x Ag = de a = Alo = & . 


gö a ədə 
y := 3t — y y Y 
y := 32 -y 
T'= —at+y ye pd 
The assignments for Ajo = $ t:=1+2%y = d i 
y:=3x—y 5 i 
“i= —rz-3y 
q:—-—z-4-3y ES 
y := 2x — y ğ i 
In general, Aa x Agpy2 = Aak46 
zo x +2y 
y:—-3r—y pida 2y 
Li Fy Li 
v= x +2y hi 2 
x := —x + 3y 
y := 3z — y 
The assignments for A4x+6 = v= —x + 3y = : 
y := 3z — y 
ET 5 
zo —x + 3y ET y 
y := 2x — y 
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The number of assignments for A4x+2 is equal to 2k + 2. 


Proposition 4. The mapping E:2" — 7", m > 2 


Fix Sle Fox + F3x3 Te... T imin 
Foxi Firo | Farş ... Em+1Tn 
E:(21,L2,T35-,8n)= | Four + Faxo + its +... + Fnit 
Fnac Fny122 + Fny203 +... Pos 1a 
such that Fa :— Fn-1 + Fn-2, Vn € Z, can be computed with m + 2 number of linear 


assignments. 


Proof. The given mapping E can be imagined as X := AX where À is a matrix such that 
whose entries satisfy the relation F, := PF, + Fn-2, Vn € Z. we construct the first two 
mappings as given below: 


£1 := £1 + 22-223 + 324 + 5z5 +... 


Lo := L1 + Lo + £3 + 274 + 3254... 


Next, we require only two elements (first and second) to compute the mapping. 


Tag: Fix + Foxo 


£4 := $3 + (Fo — Fi) zi + Fizo 


X5 := T3 T T4 


To := La T T5 


Tn t= Tn-2 + Tn-1 
Finally, we construct two mappings as follows. 
£1 = £3 — (F5 — Fi) xı + Fio 


Ta := T3 — 21 


This sequence of assignments compute mapping E effectively, and the number of these 
assignments are n + 2. 
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Lemma 5. Let E : Z? — Z? be a linear mapping defined as E(x,y) — (mx+ny, px 4- qy). 
Let E can be computed by a sequence of linear assignments 

x := aT + biy 

y := boz + aay 

go 031 + bay 


y := baz + aay 


£ := ago im + by ay 
y := byz +axy 
where ai, bi € Z, V1 «i Sk, k e Z and k > 2. Then for m, n, p, q €Z 
Mx*xQ—Px*N=M0]*02*03...0x-]1 kdk 


where m * q — pn is the determinant of the matrix ” A”, that represents coefficients of 


mapping E. 
p q 


Proof. Suppose that A denotes the determinant of the matrix A, then we have to prove that 
A = ay * a2 * 3... Qk 1 * Ak 


Consider, the evaluation of the given sequence of linear assignments 


xo a£ + biy 3.5) 

y := boat + (bəb) + a2)y 3.6) 

T :— [aza + b3(b2a1)]x a laşbı + b3(b2b1 + a2)) 3.7) 
y := [a4(b2a1) + bafaszai + b3(b2a1) lr 

+ [a4(b2b1 + az) + balaşbı + b3(b2b1 + as))Hy 3.8) 


x := [as(asa1 + b3(b301)) + bs{a4(b2a1) + balazar + b3(b201)) px 

+ las {a3b1 + b3(b2b1 + a2)) + b5{@4(b2b1 + a2) 

+ bataşbı + b3(bab1 + a2)) Huy (3.9) 
y := laş laa(b2ai) + balazas + b3(b2a1)) + befas(azar 

+ b3(bza1) + bs (a4(b2a1) + ba((asa1 + b3(b201)))]x 

+ las(aa(bəbi + a2) + ba(asbı + b3(b2b1 + a2))} 

+ be {a5(a3b1 + b3(b2b1 + a2)) + bsfar(bab1 + a2) 

+ ba (agb, + ba(bəbi + a2)) b Fly (3.10) 
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x := [ag (x-component of (k-3)th assignment) 
+ 6i, (x-component of (k-2)th assignment)]z 
+ ax-1(y-component of (k-3)th assignment) 
+ 6j, (y-component of (k-2)th assignment)]y (3.11) 
y := [az (x-component of (k-2)th assignment) 
+ b, (x-component of (k-1)th assignment)]x 
+ ax(y-component of (k-2)th assignment) 
+ by (y-component of (k-1)th assignment)]y (3.12) 


Now, 


Determinant of the Matrix — [(x-component of (k-1)th assignment) 
* (y-component of (k)th assignment)] 
— [(y-component of (k-1)th assignment) 
* (x-component of (k)th assignment)] 
—meq—pxm—A (3.13) 


Therefore, we can write 


m = la, (x-component of (k-3)th assignment) 
+ by (x-component of (k-2)th assignment)]z 
n = ax-1(y-component of (k-3)th assignment) 
+ bx 1 (y-component of (k-2)th assignment)]y 
p = [ax (x-component of (k-2)th assignment) 
+ b, (x-component of (k-1)th assignment)] = f +c 
q = az (y-component of (k-2)th assignment) 
+ by (y-component of (k-1)th assignment)] = e + d 
i (x-component of (k-2)th assignment) 


where f=a 
c = bx (x-component of (k- 


1)th assignment) 
e = aj (y-component of (k-2)th assignment) 
d = by (y-component of (k-1)th assignment) 


A=mxq—nxp=me-nf -- md -—nc 
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Since md — nc = 0, Therefore A = me — nc, that is 


A = ay. 14 [(x-component of (k-3)th assignment) 
* (y-component of (k-2)th assignment) 

+ aibi 1 [(x-component of (k-2)th assignment) 

* (y-component of (k-2)th assignment) 

— aj 184 [(x-component of (k-2)th assignment) 

* (y-component of (k-3)th assignment) 

— abı, 1 [(y-component of (k-2)th assignment) 


* (x-component of (k-2)th assignment) 
After simplification, we will have 


A = ay 184 [(x-component of (k-3)th assignment) 
* (y-component of (k-2)th assignment)] 
— [(x-component of (k-2)th assignment) 
x (y-component of (k-3)th assignment)] (3.14) 
Observe that it is again in the form of (a). Hence by continuing this process of evaluation, 


we will have 
A = ag x Gk 1 * Gk 2... 04 k aş k do * 1 


Example 16. Let E : Z? — 22 be a linear mapping defined as E(x,y) — (33x + 
307y,103x + 610y). The coefficients of mapping can be represented as a matrix A of size 


2 x 2 as follows. 
A= 33 307 
— «103 610 


Determinant (A) — —11491 
The mapping E can be computed by the following sequence of linear assignments 
:= x + 9y 
= 31 +y 
= £ + 3y 


z 
y 

z 
y:=2x+y 
c:—crcy 
y 


:= 8012x — 11491y 


Product of coefficients a;s — —11491 


51 


52 


In Situ Computation of Boolean 
Mappings 


In this chapter, we describe the in situ computation of boolean mappings. A Boolean map- 
ping describes how to determine a Boolean value output based on some logical combination 
from Boolean inputs. Boolean mappings are very important in the theory of complexity as 
well as in the design of circuits and chips for digital computers. The properties of Boolean 
mappings play a crucial role in cryptography, particularly in the design of symmetric key 
algorithms ( a class of algorithms for cryptography that use boolean function keys for both 
decryption and encryption), e.g. Two fish, Serpent, Blowfish, CAST5, RC4, TDES, and 
IDEA. 

Boolean mappings can be represented in propositional logic, or as multivariate polyno- 
mials over GF(2). To perform numerical computations on logical symbols, was one of the 
work of George Boole in the theory of logic. 

Boolean mappings have a number of applications in different areas including artificial in- 
telligence, propositional logic, circuit design, electrical engineering, game theory, reliability 
theory and combinatorics. The representation of boolean functions as a boolean polynomial 
is an extensive method in boolean algebra. Boolean polynomials play an important role 
in Reed-Muller codes (Error Correcting codes). The widespread use in electronics of inte- 
grated circuits that include “modulo 2 adders" focuses the attention on the representation 
of Boolean functions in the form of polynomials. Every boolean function can be represented 
as a polynomial. Polynomial methods have been employed extensively in circuit complexity. 
Boolean polynomials occur either directly or as a tool in the problem of decomposing a 
boolean function. The decomposition of boolean function is considered to be an important 
problem in the design of logic circuits. The decomposition of switching function is very 
important. A switching function is f : 40,1)” — (0,1). Boolean polynomial have a large 
number of applications in various fields including graph theory, law, medicine, operations 


research and spectroscopy[19, 39, 37, 3, 46]. 

Burckel et al. [11] have proved that any linear boolean mapping with input variables n can 
be computed with a double sequence of linear assignments of the same number of variables 
as given in the input. This result leads to a decomposition of boolean matrices and directed 
graphs. Burckel et al. also proved that every boolean mapping is decomposable in 4n — 3 
boolean functions. 

Next, we start by explaining the basic concepts. Then, we explain the in situ computation 
of bijective boolean mappings. Then, on the basis of this result, we explain the in situ 
computation of general boolean mappings and verified the results. We quoted definitions 
that are necessary to prove these results. At the end, we introduce a new method for the 
case of boolean bijective mappings via algebraic operations over polynomials in GF2. 


Boolean Mappings: 


We describe here the basic concept of boolean mapping and we explain it with the help of 
an example and truth table. 


Definition 4. A Boolean function/mapping is a function/mapping f from the Cartesian 
product (0, 1)" to (0,1) . Alternatively, we write f : 40,1)” — {0,1}. The set 10, 1)”, 
by definition, the set of all n-tuples (x1,--- ,æ,) where each x; is either 0 or 1, is called the 
domain of f . The set 10,1) is called the Co domain (or, sometimes, range) of f . There 
are 22”, n-ary functions for every n. 


One way to represent a boolean function whose domain is finite uses a table. Each element 
x of the domain has a row of the table listing the domain element x and the corresponding 
function value f(x). These tables make easy to understand and construct these functions 
from tables. 


Example 17. Consider the boolean mappings yı, yo, and ya from (0, 1)” to (0, 1) defined 
as follows: 


yi = Tı T T2 * T3 


Y = Lı + Lo + L1 * L3 + Li * Lo + Lo * L3 


ya = La + Tı * L2 + To * T3 


The table (4.5) represents the functions yı, yo and ya. 
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T1 | Ta | T3 vı | Y2 | Y3 
0 0 0 01010 
0 0 1 010 1 
0 1 0 0 110 
0 1 1 1 010 
1 0 0 1 1 0 
1 0 1 1 0 1 
1 1 0 1 1 1 
1 1 1 0 1 1 


Table 4.1: Truth Table 1 


4.1 Computing Bijective Boolean Mappings 


If one can compute bijective boolean mapping with an in situ program that consist of a 
set of assignments then the inverse bijective mapping can also be computed with the same 
number of assignments that can be obtained by inverting the assignments and rewriting 
from bottom to top. 


Definition 5. An in situ program that can compute a boolean mapping E : 10,1)” — 
[0,1)”, consists in a sequence of assignments of one bit component that can be written as 


tj = fj (236607) 
where f; : 40,1)” — {0,1} is a linear mapping and j is the index for the input variables. 


Theorem 1. A bijective mapping E defined over 10,1)? can be computed by an in situ 
program of the form. 


fra fee ` , İs, fas 1, 92:93:10) » In—1) Jn 
and of length 2n — 1 


Proof. We proceed by induction over n. 

For n = 1, the program will be zı := fi (zı). 

Now, we have to prove that the statement is true for n > 1. We will make use of bipartite 
graph. Suppose that G — (X, Y, A) be the bipartite multi-edges graph defined by X = Y = 
(0, 1)77İ, and (x,y) € X x Y is in A with the label zaya, for zağa € (0,1) if and only if 
E(x, £n) = (y. Yn). 

The degree of vertices of graph G will be exactly 2 due to the bijection E and graph G 
will be the union of disjoint even cycles due to the reason that it is 2-color-able regular 
bipartite graph. Therefore, by definition graph G is 2-color-able[25]. We mention that this 
is a particular case of a general result by Kónig on regular graphs [34], from which this 
proof can be generalized to any mapping on a finite set. Let us colour the edges of G with 
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elements of (0, 1). Now define two mappings E%, E! on (0,1)! and two mappings fn, gn 
from (0, 1)” to (0, 1) as follows. For each colour i € (0, 1) and every edge (x, y) with colour 
i and labeled zr, yə, define 


E'(z) — y 
Fale ta) = 
gn (y; i) = Yn 


By construction, any mapping E? is bijective on (0, 17 Then under induction hypothesis, 
each E* admits an in situ program of the form: 


enə d diues TM 
Define for every i € {0,1} and x € (0, 1)" '. 
fj (x, i) =f x) forj=n-—1,---,1 


gj (x,t) = g;(x) forj=2,---,n—1 
In other words, 


f(t, tn) = gud. a) +(1+ mie (4) forj2n-—1,--.,1 


gj(z, 24) = Ing; (x) + (1+ Ln)-9; (2) forj=n—-1,...,1 
By construction we have obtained an in situ program of length 2(n — 1)— 1 +2 = 2n — 1 
for E: 


Tn := fn(ti, mn) 
Tn-1 i= AS Tr) = Fila yib) 
Tn—2 = fy A tn) = fà xe puc) 
Ta :— fme ds = f(a, m: nea) 
T1 := fifa, yii) = Ela. ,Zn—1) 
xə — ga (1,7 En) = gis, ,En—1) 
Tra 94 0 due) 
m-ə Mo aede) 


Tn :— Qə La) = gu(yi; (0:1) 


m other words at the first step, the component x, equals a colour i, Then E'(21,--*,tn-1) 
is computed by induction in 2(n-1)-1 steps. At the step before last, we have (z1,:-- ,24 1) = 
(Y1, °°: ,Yn—1). At the last step we have x, = Yn. 
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We describe below, two examples to explain the theorem 1. 


Example 18. Consider the bijective boolean mapping E defined over (0, 11? as given below. 


21 22 T3 21 T2 T3 
0.0010 0) 0 |} 0 
0 | 0 1 1 1 1 
0 1 0 010 1 
UŞAĞI... 
1 010 1 0 1 
1 0 1 0 110 
1 1 0 1 0 | O 
1 1 1 1 1 0 


Table 4.2: Table 


The bijective mapping E1 can be computed by performing the following operations. 


a:=1+2+22%3 (fi) 
xo 1 + £1 + £2 + T1£3 + 1203 + T1T2 (F) 
13 := xı + Lo + T3 + T273 (f3) 
Lo := Li + 23 + 27123 (92) 
£1 := 1 + z1 + T2 + £3 + X33 (gi) 


After performing the first three operations fi, f}, f$ , the mapping E, transform into the 
boolean mapping that is still bijective then after performing operations g5, g, we get the 
required mapping. 


Tı | T2 | T3 Tı | T2 | T3 Tı | T2 | T3 
0 0 0 1 1 0 0 0 0 
0 0 1 1 1 1 1 1 1 
Np dı əyə ya | 00 AOE SE aa 2001 OE 
0 1 1423510 1 1 1-5”. 0 1 1 
1 0 0 1 0 1 1 0 1 
1 0 1 0 1 0 0 1 0 
1 1 0 0 0 0 1 0 0 
1 1 1 1 0 0 1 1 0 


Table 4.3: Table 


We describe another example to explain all steps in the computation of boolean bijective 
mapping. 
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Example 19. Consider the bijective boolean mapping E» defined over (0, 1)5 as given below. 


21 22 T3 21 T2 T3 
0 | 0 | 0 1 1 1 
0 | 0 1 0 1 1 
0 1 0 1 1 | 0 
ora [ects EA (E iri 
1 010 017010 
1 0 1 1 0 1 
1 1 0 1 010 
1 1 1 0 1 0 


Table 4.4: Table 


The given bijective mapping Ez can be computed by performing the following operations. 


We illustrate step by step the computation of bijective mapping E» as given below. 


= ] + £1 + T2 + £3 + 2223 


Li + Lo + T11T3 


X3 + 7172 


[= T2 23 + 1123 


:= Tı + To + L2£3 


Gü wh och 
ipa so pO dedi 


o Q 
ee Ne 


Lemma 6. Every assignment x; := f;(zxi,::: ,£n) performed in an in situ program to 
compute a bijective mapping must be linear in x;,1.e. 


fi(zi, ün En) = Ti + h(a, 19 .24—1) Tilt En) 
Proof. There are two possible cases for the mapping fi (£1, zz, £3*** , 25) for zi € {0,1} ie. 
fi (0, £2, 23 -- En) or fi (1,22,23:*- |) 


For each of these cases, there exist two possibilities for the output value i.e. 


fi (0,22,23---,24) = 0 
fi(0,z2,23:-- , £n) =1 
fi (1, 22,23:-: , £n) =0 
fi (1, £2, £3 Xa) —1 


In the computation of bijection, exactly, one of the equations from (a,b) and one of the 


equations from (c,d) will hold. 


Since the pair of equations (a,c) and (b,d) cannot hold simultaneously due to the bijection, 
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Tı | T2 | T3 Tı | T2 | T3 Tı | T2 | T3 

0 0 0 1 0 0 1 1 0 

0 0 1 0 0 1 0 0 1 

Ds Bl a an EE CR d e E E 

0 T 11510 1 11810 1 1 

1 0 0 0 0 0 0 0 0 

1 0 1 1 0 1 1 0 1 

1 1 0 1 1 0 1 0 0 

1 1 1 1 1 1 I 1 1 
T1 | T2 | T3 T1 | T2 | T3 Tı | T2 | T3 
1 1 1 1 1 1 1 1 1 
0 0 1 0 1 1 0 1 1 
“Kar | A 0: 
o eqs e la o a sı ın) 
0 0 0 0 0 0 0 0 0 
1 0 1 1 0 1 1 0 1 
1 0 0 1 0 0 1 0 0 
1 1 0 1 1 0 0 1 0 


Table 4.5: Table 


therefore, the only two possible pairs of equations (a,d) and (b,c) will hold. But the pair 
(a,d) will define the mapping 
fi (11, La, £3 s du) = Tı 
and the pair (b,c) defines 
fi (£1, La, Lg nd En) = Tı + 1 


Therefore for all cases 
fi (or, Ho , Tn) = ıı + h(x», CD Eg) 


where, h(zo,--- , 24) — 0, for the pair (a,d) 
and, h(xo,--- , 24) = 1, for the pair (b,c) 


4.2 Computing General Boolean Mappings 


In this section, we describe the in situ computation of mappings over (0,1)" with the 
supposition that two different vectors may have the same image that is not the case for 
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bijective mappings. we will describe upper bound over the length of the program that we 
use to compute these mappings. We describe some basic definitions that will help to prove 
the assertions. We will use the mapping 6 : (0, 1)” — 127) defined as 


Di, £2, T3,- cn) 2” ty +... 23 + 2.33 + a, Vn > 0 


and we also make use of simple mapping (that can be computed with n assignments ) to 
prove the upper bound over the length of program. 
The idea consist in decomposing the mapping E of (0, 1)” in Fo Po G where F and G are 
bijective and P is simple. This will lead to a sequence of 5n — 4 assignments to compute the 
mapping E. 
Definition 6. For every M > 0, a mapping f on (0,..., M} is a step mapping if for every 
0<x<y<M: 

0 x (y) - f(x) <y—x 


In particular one has f(x) € f(x +1) € f(x) +1. A mapping E on {0,1}" is simple if the 
mapping E' on [2”] such that 


E(x) = y = E'(ó(z)) = oy) 
is a step mapping. 
Proposition 5. Every simple mapping E on {0,1}", Vn > 0, is computed by an in situ 


program of length n of the form 
P1, P2,- - -Pn 


precisely, for El11,12,3,...,Un) = Vis Y2, Y3,- --;Yn and for each i = 1,2,...n 


pily, ... »Yi-1,%, E qu) — Yi 


Proof. For n input variables the n assignments are minimal. Therefore each function p; 
must return its correct final value to each of its corresponding component x;. This method 
for in situ computation of simple mapping is unique possible. However the correctness of 
this method is still remains to prove due to the fact that the mappings being computed are 
simple. 


Next, we describe a definition that will use to decompose a mapping as a composition of 
three mappings. 


Definition 7. For n > 0, let E be a mapping on 10, 1)”. A decomposition of E is a triple 
of mapping (F, P, G) on (0, 1)”, such that E = FoPoG with F, G bijective and P simple. 


For a given mapping E on ([0,1)”, the decomposition can be build in a way that we 
can group the vectors with same images via a bijective mapping G which give intermediary 
consecutive images in lexicographical order to vectors with same final images. This will 
actually maps the sets E”1(z) for x € (0, 1)" onto consecutive intervals of (27) via D o G. 
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The second step consists in the identification of vectors with the same final image via a 
simple mapping P such that 6 o P o $^! maps consecutive intervals of (27) on consecutive 
values of [2”]. The third step consists in the attribution of the correct final image value to 
each vector obtained by PoG via another bijective mapping F, that is F(P(G(x))) = Elx). 
Then F is completed arbitrarily to a bijective mapping on (0, 11^. 


Corollary 1. Every mapping E on 10,1)? is computed by an in situ program of length 
5n — 4 of the form 


f 1 1 lA 
Fi). fns Gn—1y 0665 Gly B20 ++ Pai ass 11,92) ida 


Proof. The proof can be constructed easily by using the above procedure of decomposition. 
Suppose that E = Fo PoG, i.e. the triple (F, P, G) be the decomposition of E, where 
F, G are bijective mappings and P is a simple mapping. It has been proved in theorem 1 
that both F and G can be computed by a sequence of 2n — 1 assignments respectively. 
Similarly, the simple mapping P can be computed by a sequence of n assignments as proved 
in proposition 5. Combining these number of assignments E, then can be computed by a 
sequence of 5n — 2 assignments, of the form 


, 


fise- Jn esə nəd Te reb cadis 


Two more assignments can be reduced by selecting a sequence of 2n — 1 assignments in such 
a way that it begins with the first variable for G and with the last for F. Moreover, two 
successive assignments of the same component can be combined into a single assignment. 
Therefore, gi, pı can be replaced by a single assignment gı and pn, f7, can be replaced by 
the assignment f/. 

Finally, we will get a sequence of 5n — 4 assignments. 


Definition 8. A mapping v : [2"] — N defined by v(0) +... + v(2” — 1) = 2” is called 
valuation. We denote v(l), the value of | € [2"], and by extension v(A) = Mw) for 
I€A 
A € [2"]. 
A valuation v is called boolean compatible if V i such that 0 < i < n and all 0 < j < 271, 
5 v(I) = 0 mod 2' 
j2!xl« (j-1)2* 
Definition 9. Let v be the valuation, a mapping P, : 40,1)? — 10,1)” is called projection 
of v such that for 1 from 0 to 2" — 1, v(l) consecutive elements of 10,1)? are mapped onto 
$ l(I), beginning with (0,...,0) € (0,1)" for the first | with v(l) 4 0. P, is a simple 
mapping. 
Lemma 7. Let v be a boolean-compatible valuation. 
For Li = à (127, (j + 1)2* — 1]) and for some k, k' € [271], 
where 0 € i € n and $ € [2771], we have 


Bex mp d 


k<I<k" 
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Proof. Suppose that the converse image of an interval of [2"] by the mapping ó is an interval 

of (0, 1)”. Also, the converse image of an interval (0, 1)" by the mapping P, is an interval of 

10, 1)" due to the fact that P, is a simple mapping. we have |P7 ! (1;,;)| = 5 v(l) 

j2i ie G4) 
by definition of P, 
For all | € (I,i, j), if v(I) = 0 then P; !(I;;) may be empty. we have 5 v(l) = 
jus I<(j+1)2 

0 mod 2*, due to the fact that v is boolean compatible. Thus |P*(1,¿)| = 0 mod 2. By 

induction on j for a fixed i, the result can be proved. If j = 0 then |P; ! (I; 0)| = k.2* for some 

k € [2^]. If P; 1(I;;) is not empty, then it is an interval of (0, 1)” containing (0,...,0) 

by definition of P,. Since this interval has a lenght k.2’ multiple of 27, it is of the form 
J Lu. If the property is true for all / with 0 € 1 < j, then Pl) zə) - U ra. 


0<I<k 0xi«j o<i<j’ 
Since |P, TT, yl = k.2' for some k € [277*], we must have P; !( U Li) = U lii 
O<I<j O<I<j'+k 
hence Pies) = U lii 
F<ISP ER 


Lemma 8. Let v be a boolean compatible valuation. Let a,b € 10, 1)”, and let 1 € i € m, 
If for alll > i we have a; = bi, then for alll > i we have P,(a); = P,(b)i 


Proof. We want to prove that a; = b; => P,(a)i = P, (b)i. 
Suppose that P,(a); 4 P,(b). Using lemma 7 there exist l, where, k < 1 < k’ such that 
ai ze bi. 


A contradiction. Hence our supposition is wrong. 


Proposition 6. Let E be a bijective mapping on 40,1)? computable by an in situ pro- 
gram f1,..., fn. Then the mapping E o P, is computed by an in situ program of the form 
P1,P2,.-.,Pn, Where v is a boolean compatible valuation. Precisely, for E(z1,z2,...,24) = 
(Yi, ..., gu) and for each i = 1,2,...,n: PilY1,...,Yi-1, Li, , Un) = Yi 


Proof. To prove the assertion, we use the same argument as used for the computation of 
simple mappings i.e. due to the fact that the number of assignments is minimal, necessarily 
each function p; must give its correct final value to each component x;. But the correctness 
of the method is still to be proved. 

On contrary basis, suppose that, at some step i, two different vectors z, x’ become the same 
vector, say z, whereas their final images y = E o P,(x) and y! = E o P,(z') are different. 
Because, now, they become the same vector, one has y;-1 = y; 4 = 2i-1,..., Yi = Yi = 21- 
On the other side, one has £n = zr, = Zn,...,%; = x; = zi which implies, by lemma 8 
Pula) = Pula”) = z; for all l > i. 

Assume P,(x) 4 P,(a’). Since E is bijective, then the mapping E” programmed by the 
sequence f1,...,/f;-1 is also bijective. But since E is programmed by fi,..., fn, where 
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each component is modified at most once, we must have E'(P,(x)) = E'(P,(a)) = 
(25i... 25 2_4,.--,21). A contradiction. Hence P,(x) = P,(x') and y = y 


Definition 10. A boolean compatible decomposition of E defined over (0, 1)" for n > 0 is 
a triple of mappings (F, P,G) on (0, 11^, such that E = FoPoG with F, G bijective and P 
is the projection P, of a boolean compatible valuation v = v' o c, where v'(l) is the number 
of elements of which image by E is O(1). 


We describe an algorithm in the following that build one of the mapping as described in 
the above definition. 


Algorithm: 


Define a valuation v by v(I) =|E=*(4(1))| 
Define a permutation c such that v o o is boolean-compatible 
Define a bijection G of (0, 1)” compatible with voc 
Set i=0; 
for | from 0 to 2” — 1 do 
for j from 1 to vo o(l) do 
for every x € 10, 1)" with E(x) = p(1) do 
G(x):6” (i) 
1:=1+1; 
end 
end 


end 
Define P = P,; 
Define a bijection F of (0, 1)” such that E = Fo P o G; 
for Every x € 10,1)” do 
| F(P(G())):—EQ9; 
end 
Complete the definition of F, keeping F bijective. 


Theorem 2. Every mapping E on (0, 11^ is computed by an in situ program of length 4n —3 
of the form 


one Ey dai 2000000505) 
Proof. Let (F, P, G) be a boolean-compatible decomposition of E (definition 10) . Let 
Th sg resp fp... Fm, gn ques. 91 be the sequence of 2n — 1 assignments 


computing F, resp, G, by theorem 1. Let F’ be the mapping on 10,1)” defined by the 
sequence f{,..., f/,. Then by proposition 6 the mapping F’ o P is computed by a sequence 


of n assignments f1,...,fn. Then E is computed by the sequence of 4n — 2 assignments 
fis Fe Gh a5. 9Y, fase, fas 9-1) - + -> 91, where the indices are the indices of the con- 


cerned variables. Of course gy, fı can be replaced by a single gı. Then we have a sequence 
of 4n — 3 assignments. 
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Next, we introduce a method for the case of boolean bijective mappings via algebraic 
operations over polynomials in GF(2). 


4.3 A First Tool 


In this section, we describe a method that computes a boolean function linear in zı and 
yi, where yi, Y2, ya, ..., Un are bijective mappings over £1, 19, 23, ..., En and 27, T2, 23, 

.., T, are inverse bijective mappings over yi, Ya, Y3; ..., y». In the context of Lemma 6, 
this will lead to the existence of in situ computation of bijective boolean mappings. Further 
investigation will lead to design an efficient algorithm that will help to reduce the complexity 
and ultimately to reduce the memory usage. 


4.3.1 The Algorithm 


Consider boolean bijective mappings as given below: 


yı := fi (£1, T2, 33,77 , Ty) 
Ya :— fr Gy X2; 33,777 7. 
ya = fa (£1, Ta, 23,*** ,Tn) 
UT. = nü (21, X2,X33,777 ,Tn) 
Un :— İn (21, T2, 13,777 stn) 


We are interested to compute a polynomial that is linear both in zı and yı. We can compute 
such kinds of polynomial by proceeding the following steps. 

Step-1: 

Compute the inverse mappings of each of yi, Ya, Ua. ^: ,Yn—1) Yn- These inverse mappings 
can be written as 


X1:—91 (yr, 12) V3, °°: Yn) 
X2:— 92 (yr, 12) V3, °°: Yn) 
T3 := 93 (yr, 102) 1039" Yn) 


In—1 = gai (Y1, 12) Y3,*** , Un) 
Ln := gn (Y1, Y2; Y3,"** s Yn) 
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Step-2: 
Compute the possible products of permutations of y1, Y2, ya, *** ,Un—1) Yn- 


Y1Y2,Y1Y3,*** ,Y1Yn 
Y2Y3,Y2Y4,***  U2Un 
Y3Ya4,Y3Y5),*** ,Y3Yn 


Y1Y243,Y1Y2Y4,*** ,Y1Y2Yn 
Y1Y3Y4, Y1Y3Y5)*** S V1U3Un 


Y1Y2Y3 U4, Y1Y2Y3Y5,***  V1U2U3Un 
Y1Y2Y3Y4Y5, Y1Y2Y3Y4Y6; ^: ,Y1Y2Y3Y4Yn 


Step-3: Compute the possible products of permutations of 11, £2, 23, ::: , Za 1, Tn €g. 


1172,11T3,*** ,TilTn 
123, 1284,*** , 12Tn 


Yada L385,*** ,Tazln 


111283, 11124," ** , LIT2Tn 


21834, 11 T385,*** (21730 


U1 121384), 1] L283T5,*** 2102210300 


X172732475, 1] L2L 346), ** , L1L213T4Ln 


We want to obtain a polynomial of the form 


gi + ha (Ya, yas Vasco Yn) = zi + hə (La, T3, 24,77: , En) (4.1) 
We can introduce coefficients (a1, az, a3, +*+ , da, hı, ) to generate an expression of the form 
Yı + diyə + 02Y3 + a3Y2Y3 + >: (4.2) 


Eliminate the products involving zı, and establish a system of equations in variables, 
d1,023,::: ,Q4,,7::. Substituting back these values of ar, aş, "hl ,a;,,*:- in (4.2), we will 
have the required polynomial. 
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4.3.2 Explanation and Construction 


Consider the boolean bijective mappings as given below: 


Yı := 1 + T2 + T123 


yə: 1 + £1 + Lo +21 + 2223 + $123 


Y3 := zı + Lo + T3 + $233 


These mappings can be represented in the form of a table as given below: To compute the 


Tı | T2 | T3 Yı | Y2 | Və 
0 0 0 1 1 0 
0 0 1 1 1 1 
0 1 0 010 1 
0 1 1 0 1 1 
1 0 0 1 0 1 
1 0 1 0 1 0 
1 1 0 01010 
1 1 1 1 0 0 


Table 4.6: Truth Table 2 


polynomial that is linear in zı as well as in y1, we will proceed as follows. 
Compute the possible products as given below: 


yiyə :=1+21+2%2+21%2 


Y1Y3 := Tı + T3 + zirə + 1233 


Y2Y3 := T3 T 1173 


Y1Y2Y3 :— X3 + 1273 + 1173 + 112213 
Compute the Inverse mappings: 


21 :=14+Y1Y3 + Y9Y1 + Y3 
Lo :=1+Y1Y3 + yola + Ya 
T3 := yi + Y2 + Y1Y3 


Compute the possible products: 


2112 :=1+Ya + ya + Y2Y3 

11U3 := Y1 + Y2 t Y1Y3 + Y2Y3 

T2T3 := yi + Y2Y3 + Y1Y3 + Y2Y1 
1172283 := Yı + Y1Y3 + Y2Y1 + Y1Y3Y2 
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For a, 6 and y we can make an expression 


yi + Qya + By3 + "ys (4.3) 


= (1 + X2 + 1113) 
+a (1 + 21 +22 Dir + 903 + 1123) 
+ B (£1 + £2 + 23 + 2223) 


y (za + 2123) 


We can establish a system of equations by comparing the coefficients of unwanted products 
e.g. 1173 etc. 
l+a+y=0 


a=0 
aT 8-1 
ax 


Substituting the values of a, 3 and y in (4.3), we get the required solution. Hence the 
solution that is linear in zı and y; is given as under. 


Tata + T1 := 1 + y1 yoys + ys 
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Maple Implementation 


In this chapter, we describe complete code of algorithms that we have discussed in the 
previous chapters. We implemented these algorithms in Maple language. 


5.1 Algorithm (In situ computation over Z} : 


e In this section, we describe an In Situ algorithm that can compute linear mappings 
sequentially over the set of integers Z. We split the algorithm in three parts so that it 
is convenient to illustrate. This algorithm takes a square matrix as input and return 
the corresponding sequence of assignments as out put. Moreover it return the number 
of assignments used. 


Input: Square matrix and its size 
Output: Assignments 
Assignments:—proc(); 
Local variables i, j, s; 
Global variables P, M, N; 
for i from N by —1 to 1 do 

s:=0; 

for j to N do 

| s:=s+ Mii, ils xD]; 

end 

P :— [op(P), zlil = s]; 
end 


End proc: 
Algorithm 1: First part of the Algorithm 
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Input: Square matrix and its size 
Output: Lower Triangular Matrix 
Lower Triangular:=proc(); 
Local variables i, j, a; 
Global variables P, M, N, x, X, Z; 
P:=[]; 
This section will triangulate the matrix; 
for ? from 1 to N do 

for j from i+1 to N do 

Coloperations(i,j); 
print (i,j,M); 

end 
end 
print(” Lower Triangular Matrix” ); 
print(M); 
print(” Assignments are:” ); 
Assignments(); 
print (P); 
print(” Number of assignments=” ,)nops(P)); 
print(” Verification” ); 
for ? to N do 
X [i] := Zil 
end 
for a in P do 
for i to N do 

if op(1, a) = xfi] then 

| X[i] := eval(subs(x — X, op(2, a))); 


end 

for ¿to N do 
print (X[i]); 
end 

End proc; 


Algorithm 2: Second part of the Algorithm 
We split the algorithm in three parts so that it could be adjusted on the page. Next, we 
will describe the third part of the algorithm, that take Square matrix and its size as input, 
perform column operations and return the corresponding sequence of assignments. 
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Local variables a, b, c, cc, d, dd, k, Global variables P, M, N; 
while M [i, j| 4 0 do 
a :— M{i, i|; 
b := Mii, j]; 
if b < 0 then 
P := [op(P), elil = -zlill 
for k to N do 
| Mk, j] :2 —M[k, jl; 
end 
end 
else if a < 0 then 
P := [op(P), zli] = —li]]; 
for k to N do 
| M[k, i] := —MIk, i]; 
end 
nd 
lse if a > b then 
cc := trunc(a/b); 
dd := cc — 1; 
d := a mod b; 
if d = 0 then 
P := [op P), elil = dd * ali] +11); 
for k to N do 
| Mİk,il := M[k, i] — dd x Mk, il, 
end 
nd 
lse 
P := [op(P), x[j] = cc x zlil + e[3]); 
for k to N do 
| Mik, i] := M[k, i] — cc x MIE, 5]; 


© © 


© © 


end 
lse if a — 0 then 
P := [op( P), x[] = eli] + el]l; 
for k to N do 
| Mİk,il := M[k, i] — MIE, 5]; 
end 
nd 
lse if b >= a then 
c :— trunc(b/a); P := [op(P), xfi] = x[i] + cx zİ/ll: 
for k to N do 
| Mİk, 3] :2 M[k, j] - c» Mk, i]; 


O 


© © 


end 
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5.2 Algorithm (Computing Boolean polynomial} : 


e The Algorithm construct truth tables and then bijective boolean mappings from these 
truth tables. It forms a system of equations and finally return possible polynomials 
that are linear both in zı and yi. 


N:—3 (any integer); M := 2; 
binary:— proc(a,k); 
local variables ilb; 
b:=a; 1:=[]; 
for i to k do 

l:=[b mod 2, op(1)]; 

b :— trunc(b/2); 
end 
l; 
End proc: 

Algorithm 3: A part of the Algorithm 5.2 


binaryL:= proc(a,k); 

local variables i,l,b; 

b:=a; 1:=[]; 

for i to k do 
1:=[op(1),b mod 2]; 
b:=trunc(b/2); 

end 

l; 

End proc: for ito M do 

| X[i]:=binary(i-1,N); 
end 


End proc: 
Algorithm 4: A part of the Algorithm 5.2 
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5.2. ALGORITHM {COMPUTING BOOLEAN POLYNOMIAL} : 


sous:=proc(f, p); 
Local variables a, ok; 
ok:=false; for a in f do 
if convert(expand(a + p), mod2) = 0 then 
| ok:=true; 
end 
end 
ok; 


End Proc: 
Algorithm 5: A part of the Algorithm 5.2 


fctx:= proc(1); 
local variables i, j, k, s, p; 
global variable X; 


s:=0; 
for ? fo M do 
if l[i] = 1 then 
p:=1; 


for k from 1 to N do 
if Xfi]k] -1 then 
p:=pxxl[k] else 
| pi- ps (1+ a[k]) 
end 
end 
end 
S'=S+p+p*xs;: 
end 


end 
s:=convert(expand(s), mod2); 
S, 
End proc: 
Algorithm 6: A part of the Algorithm 5.2 
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fcty:=proc(l): 

Local variables i, j, k, s, p; 

Global variable Y; 

s:=0; 

for ? fo M do 

if {/ij=1 then 

p:=1; 

for k from 1 to N do 

if Y/i//k/=1 then 
p := p * y[k]; 
else 

| p:=px*(1+ylkl) 

end 


end 

end 
s:=sSs+p+p*s; 
end 


end 
s := convert(expand(s), mod2); 
s; 
End proc: 

Algorithm 7: A part of the Algorithm 5.2 
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test:=proc(f) ; 

for i to M do 

| Ylij:2binary (ffi],N); 
end 


print (^ 
print(f); 
for i to M do 
print(X [i], Y [1]); 
end 
for j to N do 
L-[; 
for i to M do 
| k-l[ep(D, YE: 
end 
YY [j]:=fctx(1); 
end 
for j to N do 
L=||; 
for ? fo M do 
| 1—lop() Xİ 
end 
XX[jl:-fety(1): 
end 
print(” MAPPING"); 
s:=[]; 
for j to N do 
| s:=[op(s),y [il]; 
end 
ss:={}; 
for j to N do 
ss:=ss union (y (¡]=YY[j]); 
end 
for k to N do 
if k=2 then 
| print(” products"); 
end 
t:=combinat [choose] (s,k); 
for pp in t do 


p:=1; 

for a in pp do 
| P:=3Pp*a; 

end 


q:=subs(ss,p); 
q:=convert(expand (q),mod2); 
print(p,”===”,q); 

end 


end 


Algorithm 8: A part of the Algorithm 5.2 
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print(” INVERSE MAPPING"); 
s:=l); 
for j to N do 

| s:=[op(s) xi] 
end 
ss:={}; 
for j to N do 
ss:=ss union (x|¡]=XX[j]); 
end 
for k to N do 
if k=2 then 

| print(” products"); 
end 
t:=combinat [choose] (s,k); 
for pp in t do 


p:=1; 

for a in pp do 
| P:=3Pp*a; 

end 


q:=subs(ss,p); 

q:=convert (expand(q),mod2); 
print(p," ===” ,q); 

end 


end 

print(” SOLUTIONS THAT ARE LINEAR IN Y1 and X1"); 
col([0]); print” SOLUTION BY LINEAR SYSTEM"); 
syst(): 


End Proc: 
Algorithm 9: A part of the Algorithm 10 


prx:— proc(k): 
Local variables i,l,p; 
Global variables N; 1:-binaryL(k,N): 
p:=1; 
for i to N do 

if l[i] = 1 then 

| p:=px*zli); 
end 


Algorithm 10: A part of the Algorithm 5.2 
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syst:=proc(); 

Local variables Z, i, 1, p, px, py, pyx, 5, c, k; 
Global variables XX,Y Y,M; 

s:— (kiyfl]] = 1): 

print(” CONDITIONS"); 

print(” MUST BE LINEAR IN Y1==¿”,s); 
S:— 8; 

for i to M-1 do 

px:=prx(i); 

if has(px,x/1]) then 

c:=0; 

for j to M-1 do 

py:=pry (j); 

pyx:=pryx(j); 


| ecc k[py]; 
end 
end 
if pr = z[1] then 
c:= (c= 1); 
else 


end 
if c z (0 = 0) then 
8 := sunionc; 
end 
if pz = 2/1] then 
cond:=” MUST CONTAIN” ; 
else 

| cond:="MUST HAVE NO” 
end 


? 
end 
print(cond, px,” ===j”,c); 
end 
end 


Algorithm 11: A part of the Algorithm 5.2 


if sous(pyz + Z,px)and((py € y[1])ornothas(py, y[1])) then 
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print(” THE LINEAR SYSTEM IS”); 

print(s); 

s:=msolve(s,2); 

print(” THE GENERAL SOLUTION FORM IS”); 


print(s); 
s := subs(z1 = 0, s); 
s := subs( z2 = 0, s); 
s := subs( z3 = 0, s); 
s := subs(z4 = 0, s); 
s := subs( z5 = 0, s); 
s := subs( z6 = 0, s); 
s := subs( z7 = 0, s); 
s := subs( z8 = 0, s); 
print(” A PARTICULAR SOLUTION IS”); 
print(s); 
p:=0; 
for j to M-1 do 
py:=pry (j); 


c:=subs(s,k[py]); 
if cz klpy] then 
| p:= p+ cx py; 
end 
end 
print(” SOLUTION FOR THE POLYNOMIAL LINEAR IN XI and Y1=”,p); 


End Proc: 
Algorithm 12: A part of the Algorithm 13 


idem:=proc(p,q) ; 

Local variables i,ok; 

ok:=true; 

for i from 2 to N do 
if (pi Z qfil) then 
| ok:=false; 
end 

end 

ok; 


End Proc: 
Algorithm 13: A part of the Algorithm 5.2 
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col:=proc(c) : 

Local variables i, j, k, ok; 
k:=nops(c); 

ok:=true; 

for ¿to k do 

for j from i41 to k do 


if idem(X[i], X[j]) and (c[i]-c[j]) then 


| ok:=false; 
end 
if idem(Y[i], Y[j]) and (c[i]-c[j]) then 
| ok:—false; 
end 
end 
end 


if not ok then 
else if k=M then 
print(fctx(c)," === ,fcty(c)); 
else 
col(fop(c),0]): 
col([op(c).1]): 
end 
end 


end 


End Proc: 
Algorithm 14: A part of the Algorithm 5.2 


pry:=proc(k); 
Local variables i, 1, p; 
Global variable N; 
l:=binaryL(k,N); p:=1; 
for ¿to N do 

if {/ij=1 then 

| p:=px*yll; 

end 
end 
p:=convert(expand (p),mod2); 
P; 
End proc 

Algorithm 15: A part of the Algorithm 5.2 
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pryx:=proc(k); 

Local variables i, l, p; 

Global variables N, YY; l:-binaryL (k,N); p:=1; for i to N do 
if ////-1 then 
| p:=p*YY [i]; 
end 

end 

p:—convert(expand(p),mod2); p; 


End Proc: 
Algorithm 16: A part of the Algorithm 5.2 


all:=proc(f) ; 
Local variables i; 
Global variables M; 
if nops(f)-M then 
test(f) else 
for i from 0 to M-1 do 
if not member(i,f) then 
| all(İop(f).il): 
end 
end 
end 


end 


End Proc: 
Algorithm 17: A part of the Algorithm 5.2 


5.2. ALGORITHM {COMPUTING BOOLEAN POLYNOMIAL} : 


ex:=proc(k) : 

Local variables fa; 

Global variables M; 

if k=0 then 

all((D: 

if k=1 then 

f:=[]; 

while nops(f) < M do 
a:=rand() mod M; 
if not member(a,f) then 

| f:=[op(f),al; 

end 

end 

test(f); 

end 


end 


End Proc: 
Algorithm 18: A part of the Algorithm 5.2 
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5.3 Algorithm {In situ computation over Z/NZ} 


e Algorithm that compute a mapping E over Z/NZ by a sequence of linear assignments. 
This algorithm takes a square matrix as input and return the corresponding assignment 
matrices as out put. We can construct sequence of linear assignments directly from 
these assignment matrices. 


gcd-multipliers := proc(x::list, alpha::integer, i::integer); 
Local variables tmp , fac, projections,1,j,k,mres; 
fac := ifactors(alpha) [2]; 
projections : []:; 
for k from 1 to nops(fac) do 
l := [seq(0, j—1..nops(x))]; 
lfi] := 1; 
if igcd(x[i], fac|k]1]) 4 1 then 

for j from 1 to nops(x) do 

if igcd(z[j], fac[k]1]) = 1 then 
il = 1; 


15) := 1; 
break; 
end 
end 
end 
projections := [op(projections), 1); 
end 
tmp :— [seq(fac[k,1]^fac[k,2], k=1..nops(fac))]; 
mres : İl: 
mres :=[op(mres), chrem(projections, tmp)][]; 
End Proc: 


Algorithm 19: A part of the Algorithm 5.3 


creatematrix :— proc(d::integer); 
Local variables M, i; 
M := Matrix(d,d,0); 
for i from 1 to d do 
| Mfi,i] := 1: 
end 
return M; 


End Proc: 
Algorithm 22: A part of the Algorithm 5.3 
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5.8. ALGORITHM (IN SITU COMPUTATION OVER Z/NZ} 


isone :— proc(M::Matrix); 
Local variables i, j; 
for i from 1 to RowDimension(M) do 
for j from 1 to ColumnDimension(M) do 
if Mİ, j]! = 1 then 
| return false: 


end 
return true; 


End Proc: 
Algorithm 20: A part of the Algorithm 5.3 
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with(List Tools); 

matrixop := proc(Ma::Matrix, N::integer); 

Local variables i, j, k, nr, nc, left, right, g, 1, T, G, M, U, final, F; 
M := copy(Ma); 

left := []; 

right := []; 

nr := RowDimension(M); 

nc := ColumnDimension(M); 


for k from 1 to min(nr, nc) do 
g:— iged ( seq(Mİi,kl, i—1..nr) ); 


if g — 1 then 
continue: 
end 


if igcd( Mik, kl/g, N) 4 1 then 

1:— gcd-multipliers([seq(M [i,k]/g, i=1..nr)], N, k); 
T := creatematrix(nr)mod N; 

for j from 1 to nr do 

| Tli — Hi: 

end 

M := Multiply(T , M)mod N: left :— [op(le ft), Matriz Tnverse(T)modNl, 
end 

G:= creatematrix(nc); 

G[k,k] := M[k,k]: M := Multiply(M , MatrixInverse(G))mod N; 
U := creatematrix(nc); 

for j from 1 to nc do 

| Ulk, il: Mİk, j]; 

end 

M := Multiply(M , MatrixInverse(U))mod N; 

right := [op(right), Multiply(U , G)mod N |; 

end 

final := |]; 

for i from 1 to nops(left) do 

if isone(left/i]) = false then 

| final :— [op(final), left [i] 1: 

end 


end 

for i from nops(right) to 1 by -1 do 
if isone(right/i]) = false then 

| final :=[op(final), rightlill: 
end 


end 
Reverse(final); 


End Proc: 
Algorithm 21: A part of the Algorithm 5.3 


83 


5.4. ALGORITHM: 


5.4 Algorithm: 


e Algorithm that compute a mapping E : Z" —> Z", m > 2, where Z is a set of 


integers, defined as 


E:(xi,22,93,:7 ün) m 


such that Fu := 


Fix, + Foto + Far +... Fas 

Fox, Firo Farş Fe se Fm+1Tn 

Pari + Fazo + F5£3 +... + FEm+2Tn 
Ent1 + Fny122 + Fn+283 +... Posi 


n—1 + Fn-2, Vn € Z, can be computed with m + 2 number of linear 


assignments. The algorithm takes only 2 initial entries of the square matrix A at positions 
aıı and dış and construct all other coefficients and return the linear assignments that 


can compute the mapping E. 


Input: integer 
Output: fibonacci numbers 
oo:— proc(n); 
local variables i, f, s; 
s:=[1); 
ffl] := 1;f[2] = 1; 
for i from 2 to n do 
fli + 1) = fli + fl 1); 
s :— [op(s), fli + 11); 
end 
End proc; 


Algorithm 23: a part of the Algorithm 
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Input: integer 
Output: assignments 
ass:= proc(N): 
local variables j, jj, ii, xxx, first, second, m, p, PP, S, xx; 
p:—0; pp:=0; first:=2; second:=3 ; 
if N < 2 then 
print(”N must be greater than 2”): 
else 
s:=00(N); 
for j from 2 to N do 
| p:=p+s[j-1]*x[j]; 
end 
xx[1]:=x[1]+p; 
for jj from 2 to N-1 do 
| pp:=pp+s[ij-1)*x[jj+1); 


end 

print(x[1]," =” xx[1]); 

xx[2]:=xx[1]+x[2]+pp; 

print(x[2]," =” ,x[1]--x[2]--pp) 

xx[3]:=first*xx[1] +second*xx[2]; 

print(x[3]," =” first*x[1] +second*x[2]); 

if N > 3 then 
xx[4]:2xx[3]-- ((second-first) *xx[1] +first*xx[2]); 
print (x[4],” =” ,x[3]+ (second-first) *x[1]--first*x[2]); 


if N > 4 then 
for iii from 5 to N do 
| xx[iii]:=xx[iii-2]4-xx[iii-1]; 


end 

end 
end 
for v from 5 to N do 

print (x[v],” =” x[v-2]+x[v-1]); 
end 
print(x[1]," :=” ,x[3]-((second-first) *(x[1])+(first) *(x[2]))); 
print(x[2]," :—" ,xİ31-xİ11): print (? Pe eee aaa AAA A ) 


xxx[1]:=xx[3]-((second-first)*(xx[1])+ (first) *(xx[2])); xxx[2]:=xx[3]-xxx[1]; 
print(x[1]=xxx[1]); print (x[2]=xxx[2]); 

for ii from 3 to N do 

| print (x[ii]=xx{ii]); 


end 
end 
end 
End; 
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5.5 Conclusions: 


The in situ computation of mapping by a sequence of assignments, without using any 
extra variable other than the variables available as input variables, is a new way in code 
optimization according to the best of our knowledge. This way of computation helps 
directly to reduce the memory usage and ultimately the idea can be implemented in 
compiler / processor optimization, register allocation and chip designing. 

The sequential decomposition of matrices is a powerful tool for computations of linear 
transformations using minimal memory. We prove that, for such computation extra 
memory than input variables is not required. We prove the sequential decomposition of 
linear mapping over field and extended this idea over the Ring and set of integers. It is 
still an open path for the researchers to improve the efficiency of these algorithms and to 
find the minimum number of assignments required to compute a linear mapping using in 
situ computation. The parallel applications of this work, in the decomposition of matrices 
may lead to solve system of linear equations, to triangularize a matrix efficiently etc. 
Since the matrices can be triangulate with different methods, therefore we can find the 
applications of in situ computations of mappings corresponding to these methods. For 
example, Triangularization of a square matrix, using Extended Euclidean Algorithm leads 
to find the sequence of linear assignments that can compute linear mappings using in situ 
program, consists of these linear assignments. The upper bound, of Extended Euclidean 
Algorithm, that has been found by Lame theorem can be helpful to find the number of 
assignments used by corresponding in situ program, and the lower bound will yield a new 
interesting results both for such algorithms and in situ programs. 


We have also described that in situ computation of general and bijective boolean map- 
pings and presented a suitable bound for this computation. We introduce a new method on 
the basis of bijective boolean mappings via algebraic operations over polynomials in G F2. 
This new method is a powerful tool to compute in situ boolean mappings. Initially, this 
method compute linear boolean polynomials that are linear both in zı and yi. 

'The main motivation of the work presented is the processor performances because the ap- 
plication of this work concerns to hardware as well as software in the sense of chip designing 
and compiler optimization. 
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